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Abstract 

Using the matrix Ricmann-Hilbert factorization approach for nonlinear evolution sys- 
tems which take the form of Lax-pair isospectral deformations and whose corresponding 
Lax operators contain both discrete and continuous spectra, the leading-order asymp- 
totics as t —^ ±oo of the solution to the Cauchy problem for the modified nonlin- 
ear Schrodinger equation, idtu+ ^dlu+\u\'^u + isdx{\u\'^u) — 0, s £ M>o, which is 
a model for nonlinear pulse propagation in optical fibers in the subpicosecond time 
scale, are obtained: also derived are analogous results for two gauge-equivalent nonlin- 
ear evolution equations; in particular, the derivative nonlinear Schrodinger equation, 
idtq+d^q—idx{\q\^q) = 0- As an application of these asymptotic results, explicit expres- 
sions for position and phase shifts of solitons in the presence of the continuous spectrum 
are calculated. 
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1 Introduction 



With the current emphasis on the utiUzation of optical fibers, capable of supporting soli- 
tons, as the communication channel in the practical realization and implementation of 
all-optical (lightwave), ultrahigh-bit-rate, long-distance communication systems using the 
return-to-zero (RZ) format for generating the optical bit stream, design issues requiring 
the consideration of several factors, e.g., soliton widths and inter-soliton spacings, are 
intimately related to the study of the fundamental dynamical processes associated with 
the propagation of high-power ultrashort pulses in optical fibers (at the present stage of 
technology, these systems can at best still only be called near-soliton(ic)-based) Q. The 
standard, classical mathematical model for nonlinear pulse propagation in the picosecond 
time scale in the anomalous dispersion regime in an isotropic, homogeneous, lossless, non- 
amplifying, polarization-preserving single-mode optical fiber is the nonlinear Schrodinger 
equation (NLSE) [||, |2|. However, experiments and theories on the propagation of high- 
power ultrashort pulses in the subpicosecond-femtosecond time scale in monomode optical 
fibers have shown that the NLSE is no longer a valid model, and that additional nonlinear 
terms (dispersive and dissipative) and higher-order linear dispersion must be taken into 
account: in this case, pulse-like propagation is described (in dimensionless and normalized 
form) by the following nonlinear evolution equation (NLEE) 

1 ~ o r 

idzU H — d'^u + \u\'^u + isdri\u\'^u) = —iVu -\- i5d^u + -^udr{\u\^), (1) 

2 To 

where u is the slowly- varying amplitude of the complex field envelope, z is the propagation 
distance along the fiber length, r is the time measured in a frame of reference moving with 
the pulse at the group velocity, s G M>o governs the effects due to the intensity dependence 
of the group velocity (self-steepening), F is the intrinsic fiber loss, 5 governs the effects of 
the third-order linear dispersion, and t„/to governs the soliton self-frequency shift effect. 

Since, under typical operating conditions, F, ^, and tVi/tq are small parameters 0], a 
strategy to study the solutions of Eq. (1), for which the nonlinear effects dominate the 
higher-order linear dispersive one, is to set the right-hand side equal to zero, thus obtaining 
the following NLEE (integrable in the sense of the inverse scattering method (ISM) ) , 

idtu + ^d'^u + \ufu + isdx{\u\'^u) = 0, (2) 

which, hereafter, is called the modified nonlinear Schrodinger equation (MNLSE) (the phys- 
ical variables, z and r, have been mapped isomorphically onto the mathematical t and x 
variables, which are standard in the ISM context), and to treat Eq. (1) as a non-integrable 
perturbation of the MNLSE. From the above discussion, it is clear that perturbations of 
multi-soliton solutions of the MNLSE can be very important in the physical context related 
to optical fibers Q. Since practical lasers excite not only the soliton(ic) mode(s), but also an 
entire continuum of linear-like dispersive (radiative) waves, to have physically meaningful 
and practically representative results, it is necessary to investigate solutions of the MNLSE 
under general initial (launching, in the optical fiber literature |]l|) conditions, without any 
artificial restrictions and/or constraints, which have both soliton(ic) and non-soliton(ic) 
(continuum) components: it is towards such a solution that the initial pulse launched into 
an optical fiber is evolving asymptotically 0]. In physical terms, the pulse adjusts its width 
as it propagates along the optical fiber to evolve into a (multi-) soliton, and a part of the 
pulse energy is shed away in the form of dispersive waves in the process: normally, these 
dispersive waves form a low-level broadband background radiation that accompanies the 
(multi-) soliton 0, |^. From the physical point of view, therefore, it is seminal to understand 
how the continuum and the (near-) soliton(s) interact, and to be able to derive an explicit 
functional form for this process. Since Eq. (2) is integrable via the ISM, one can use one 
of the techniques developed in the framework of this approach to solve the aforementioned 
problem; in particular, in this paper, the Riemann-Hilbert (RH) factorization method is 
applied. 
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For several soliton-bearing equations, e.g., KdV, Landau-Lifshitz, NLS, sine-Gordon 
and MKdV, it is known that the dominant (0(1)) asymptotic (t — > ±00) effect of the 
continuous spectra on the multi-soliton solutions is a shift in phase and position of their 
constituent solitons [||: as will be shown in this paper, an analogous, though analytically 
more complicated, situation takes place for the MNLSE (the additional complexification 
occurs due to the non-standard normalization of the associated RH problem). While the 
above-mentioned works deal only with the leading-order (0(1)) asymptotic term, in this 
paper, for the MNLSE, not only the leading-order, but the next-to-leading-order (0(t~^''^)) 
term as well is derived; in particular, besides inducing an 0(1) position and phase shift on 
the multi-soliton solution, this 0(t~^/^) term represents the evolution of the continuum 
component (dispersive wavetrain [|l|) as well as the non-trivial interaction (overlap) of the 
soliton and continuum components of the solution. It is worth mentioning that, even though 
there have been several papers Q devoted to studying the soliton solutions of the MNLSE, 
to the best of our knowledge, very little, if anything, was known about its solution(s) for the 
class of non-reflectionless initial data until very recently |^. In the framework of the ISM, 
an asymptotic analysis of the aforementioned solution for the MNLSE can be divided into 
two stages: (1) the investigation of the continuum (solitonless) component of the solution 
P, |9|, and (2) the inclusion of the (multi-) soliton component via the application of 



a "dressing" procedure |11, 12 1 to the continuum background. In this paper, the above- 
mentioned asymptotic paradigm is carried out systematically for the MNLSE: the results 
obtained in this paper are formulated as Theorems 2.1-2.3. 

This paper is organized as follows. In Sec. 2, a matrix RH problem for the solution of a 
NLEE gauge-equivalent to Eq. (2) is stated, and the results of this paper are summarized 
as Theorems 2.1-2.3. In Sec. 3, an extended RH problem is formulated and shown to be 
equivalent to the original one stated in Sec. 2, and as it is shown that the solution 

of the extended RH problem converges, modulo exponentially decreasing terms, to the 



solution of a model RH problem. In Sec. 4, the Beals-Coifman [13| formulation for the 
solution of a RH problem on an oriented contour is succinctly recapitulated, and the model 
RH problem is solved asymptotically as t ^ -|-co for the Schwartz class of non-reflectionless 
generic potentials. In Sec. 5, a phase integral which is associated with the non-standard 
normalization of the above-mentioned RH problem is evaluated asymptotically as t^-\-oo. 
Finally, in Sec. 6, the asymptotic analysis as t— >— 00 is presented. 

2 The Riemann-Hilbert (RH) Problem and Summary of 
Results 

In this section, the matrix RH problem is stated, and the main results of the paper are 
formulated as Theorems 2.1-2.3. Before doing so, however, it is necessary to introduce 
some notation and definitions which are used throughout the paper: 

Notational Conventions 

(1) ^a(3, G {1;2}, denote 2x2 matrices with entry 1 in (a/3), {ea(3)ij '■= ^ai^pj, 
i,jG{l,2}, where 5ij is the Kronecker delta; 

(2) I: = eii-|-e22 = diag(l, 1) denotes the 2x2 identity matrix; 

(3) C73:=eii-e22 = diag(l,-l), C7_:=e2i, cr+:=ei2, and ai:=a^+a+; 

(4) for a scalar vo and a 2x2 matrix T, zu^^^'^^^T ■.=tu'^^Tzu^'^^; 

(5) (•) denotes complex conjugation of {•); 

(6) M2(C) denotes the 2x2 complex matrix algebra with the following inner product 
((•,•): M2(C)xM2(C)^C), Va,6£ M2(C) , (a, &) := tr(&a), and (foraeM2{C)) the 
norm on M2{C) is defined as \a\ := ^ {a,a); 
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(1) CPiD;M2{C)) := {/; f:D^ M^iC), \\f\\c,^D-M,iC)) ■= UDlfioriM)'^" < oc,p G 
{1,2}}; 

(9) for D an unbounded domain o/ M U zM, let S{D;C) (resp. S{D;M2{C))) denote the 
Schwartz class on D, i.e., the class of smooth C-valued (resp. M2(C) -valued) functions 
f{x):D^C (resp. F{x): D ^ M2{C)) which together with all derivatives tend to zero 
faster than any positive power of \x\~^ as |x| — >co. 

In this paper, as in ||^], along with the MNLSE, the following NLEEs are studied: 

idtQ + d^Q + iQ^d^Q + lQ\Q\^ = 0, (3) 

with initial condition Q{x,0) g5(M;C), and the derivative nonlinear Schrodinger equation 
(DNLSE), 

idtq + d^,q-idM^q)=0, (4) 

with initial condition q{x,0) £ 5(M;C). To recall the relations between the solutions of these 
NLEEs, the following propositions are formulated: 

Proposition 2.1 (|l^) The necessary and sufficient condition for the compatibility of the 
following system of linear ODEs (the Lax pair) for arbitrary AgC, 

t;A) = U{x, t;A)^'(x, t;A), t;A) = V{x, i;A)^'(x, t;A), (5) 

where 

U{x, t;X) = -iAV3+A(ga_ +ga+) - flQI Vg, 
y(x,t;A) =2A2[7(x,t;A)-a((a,g)a„-(9,.Q)a+) + (i|Q|4 + i(g5,.Q-Q9,.Q))a3, 
is that Q{x,t) satisfies Eq. (3). 

Proof. Eq. (3) is the Frobenius compatibility condition for system (5). ■ 

Proposition 2.2 Let Q{x,t) be a solution of Eq. (3). Then there exists a corresponding 
solution of system (5) such that ^'(x,t;0) is a diagonal matrix. 

Proof. For given Q{x,t), let ^{x,t;X) be a solution of system (5) which exists in ac- 
cordance with Proposition 2.1. Setting A = in system (5), one gets that ^{x,t;X) = 
exp{—^J^JQ[g,t)\'^dg}}C{X), for some xq € R and non-degenerate matrix IC{X) which is 

independent of x and t. The function ^(x,t;A) := ^(x, t;A)(A^(A))~^ is the solution of sys- 
tem (5) which is diagonal at A = 0. ■ 



Proposition 2.3 (|T^) Let Q{x,t) be a solution of Eq. (3) and ^{x,t;X) the corresponding 
solution of system (5) given in Proposition 2.2. Set '^q{x,t;X) :='if~^{x,t;0)'if{x,t;X). Then 

d^^q{x, t;X) = Uq{x, t;A)^'g(x, t;X), 9t^g(x, t;A) = Vg(x, t;X)'i'g{x, t;X), (6) 

where 

Ug{x, t;X) = -iX'^as + X{qa- + qa+), (7) 

/ -2iX^ - iX^\q\^ 2X^q + iXd,q + X\q\^q\ . 
)/,ix,t,Aj- 1^ 2A3g-a9.g + A|g|2g 2iX' + iX^\q\^ J' 

with q{x, t) defined by 

q{x,t) :=Q(x,^)((^-l(x,^;0))n)^ (9) 
is the "Kaup-Newell" |16| Lax pair for the DNLSE. 
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Proof. Differentiating *q(x,t;A) := ^~^{x,t;0)"^{x,t;X) with respect to x and t and 
using the fact that ^{x,t;0) = exp{-^J^JQ{Q,t)\'^dQ}, for some xq G M, and ^{x,t;X) 
satisfy system (5) for A = and AgC\{0}, respectively, defining q{x,t) as in Eq. (9), one 
gets that ^q{x, t;X) satisfies system (6), where Uq{x, t;X) and Vq{x, t;X) are given by Eqs. (7) 
and (8): Eq. (4) is the Frobenius compatibihty condition for system (6). ■ 

Proposition 2.4 Ifq{x,t) is a solution of the DNLSE such that g(x,0) G(S(R;C), then 

<x, t) := ^ exp{i(x - - x, |) (10) 

satisfies the MNLSE with initial condition 'u(x,0) G(S(]R;C). 

Proof. Direct substitution. ■ 

Remark 2.1 A convention is now adopted which is adhered to sensus strictu throughout 
the paper: for each segment of an oriented contour, according to the given orientation, the 
"+" side is to the left and the "— " side is to the right as one traverses the contour in the 
direction of the orientation; hence, (•)+ and (•)_ denote, respectively, the non-tangential 
limits (boundary values) of (•) on an oriented contour from the "+" (left) and "— " (right) 
sides. 

Before stating the matrix RH problem which is investigated asymptotically (as t— >zboo) 
in this paper (see Lemma 2.1), it will be convenient to introduce the following notation: 

let Zd := u£i({±Ai}U{±Al}) and f := {A; '^{X'^) = 0} (oriented as in Fig. 1) denote, 
respectively, the discrete and continuous spectra of the operator dx — U{x,t;X), and cr£ := 
Spec{dx-U) = ZdUf (Z^nf = 0). 





0(A) 


3(A2)<0 


3(A2)>0 




+ SR(A) 


+ 






■~ 0(A2)<0 



Figure 1: Continuous spectrum T. 

Lemma 2.1 Let Q{x,t), as a function o/x, g5(M;C). Set m{x,t;X) :='^{x,t;X) cx];){i{X'^x+ 
2A*^t)(j3}. Then: (1) the bounded discrete set is finite (card(^d) <oo); (2) the poles of 
m{x,t;X) are simple; (3) the first (resp. second) column of m{x,t;X) has poles at {^Xi}^^ 
(resp. {zbAll^i); <^^d (4) VteR the function m(x, t;A): C\ (Z^uf ) ^ SL(2,C) solves the 
following RH problem: 

a. m{x,t;X) is meromorphic ^ XE:C\T; 

b. m{x, t;X) satisfies the following jump conditions, 

m+(x,t;A) = m-{x,t;X)v{x,t;X), AgF, 

where 

v{x, t;X) = exp{-z(A2a;+2A4t)ad(a3)} ( ^ ~ ''i^^'^^ ""^^^ ] , 

r(A), the reflection coefficient associated with the direct scattering problem for the 
operator dx — U{x,t]X), satisfies r(—X) = —r(X), and r(X) &S(r;C); 
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c. for the simple poles of m{x,t;X) at {±Xj}jLi and {±Xj}jLi, there exist nilpotent 
matrices {vj{x,t)a^}^^i and {vj{x,t)aj^}^^-^, respectively, such that the residues, for 
satisfy the (Beals-Coifman H^) polar conditions. 



res(m(x, f;A);Aj) = lim m(x, t;A)T;j(x, t)(T_, 

res(m(x, t;A);— Aj) = — (T3res(m(x, t;A);Aj)cr3, 
res(m(j;, t;A);Aj ) = lim_?7i(x, t;A)vj(x, t)(T+, 

res(m(x, t;A);— Aj ) = — cr3res(m(x, t;A);Aj )(T3, 

where t) :=Cj exp{2i(Aja;+2Ajt)}, andCj are complex constants associated with 
the direct scattering problem for the operator dx — U{x,t;X); 

d. as X^oo, X£C\{ZdlJV), 

m{x,t;X) =1 + 0{X-^). 

Proof. Conditions (l)-(4) follow from the results given in 17]. ■ 

Lemma 2.2 Let ll'^ll^oo^]^.^-) ^^PasM l^^'^)! ^ ^' ^^s'^- (V the RH problem formulated 
in Lemma 2.1 is uniquely solvable; (2) t;A) = m(x, i;A) exp{— i(A^a; + 2A^i)cj3} is the 
solution of system ( 7) with 

Q{x, t) := 2i lim {Xm{x, t;X))i2; (11) 

(3) the function Q{x,t) defined by Eq. (11) satisfies Eq. (3), and 

q{x, t) := Q{x, t){{m-'^{x, t; 0))nf (12) 
satisfies the DNLSE; and (4) m(x,t;X) possesses the following symmetry reductions, m{x,t; 
X) = a3m{x,t;—X)a3 and m{x,t]X)=aim{x,t;X)ai. 

Ifr{X)£S{T;C), then, for any t gM, Q{x,t) (resp. q{x,t)), as a function of x ,£ S {M.;C) . 

Proof. The solvability of the RH problem (formulated in Lemma 2.1) is a consequence 



of Theorem 9.3 in ]|l^ and the vanishing winding number of 1 — r(A)r(A), /j;d(arg(l — 
r(A)^)) = E;e{>}s(On(0=0, where s(>) = -s(<) = 1, and n(>) :=card({Aj; 9(A|)>0}), 

items (2) and (4) can be verified through straightforward calculations, and the fact that 
q{x,t) (Eq. (12)) satisfies the DNLSE follows from Proposition 2.3 and the definition of 
m{x,t]X). ■ 

Remark 2.2 In fact, in this paper, the solvability of the RH problem for 1 1'^l l£oo(]^.C) <1 
is proved for all sufficiently large \t\: the solvability of the RH problem for ]]r]] < 1 
in the solitonless sector, Zd = 9, for all sufficiently large \t\ was proved in Note : the 
condition ]]?']]£oo(f) ^ ^ which appears in Q is restrictive, and can be replaced by the weaker 
condition Ilr-]]£^(K.Q<1. 

Before summarizing the main results of this paper, namely. Theorems 2.1-2.3, some 
further preamble is required: (1) the Kaup-Newell parametrization jl6] is adopted for the 
discrete eigenvalues, Aj := Aj exp{|(7r— 7^)}, Aj>0, 7jG(0,7r), l<j<N, and Xj ■.=^j+ir]j, 

where = — A|cos7j and r]j = A|sin7j (note that, with this parametrization, {±Xi}fLi 

(resp. {±Xi}^Li) lie in the 1st and 3rd quadrants (resp. 2nd and 4th quadrants) of the 
complex plane of the auxiliary spectral parameter, A); and (2) it is supposed throughout that 
Ci¥'Cj! ^ ¥^ j ^ N J so that it is convenient to choose the following enumeration for the 
points of the discrete spectrum ( ordering of the solitons), £,i> ■ ■ ■ > £,n> ' ' ' > ■ 

Remark 2.3 Even though the "symbol" ("notation") C{z) appearing in the various final 
error estimations is not the same and should properly be denoted as Ci{z), 6*2(2;), etc., the 



6 



simplified "notation" C{z) is retained throughout since the principal concern here is not its 
explicit functional z-dependence, but rather, the functional class(es) to which it belongs. 
Throughout the paper, Mg]R>o is a fixed constant. 

Remark 2.4 In Theorems 2.1-2.3 (see below), one should keep the upper signs as t^+oo 
and the lower signs as t^—oo everywhere. 

Theorem 2.1 Let m{x,t;X) be the solution of the RH problem formulated in Lemma 2.1 
with the condition I kl l£oo(]^.iJ^) <1 o,nd Q{x,t) be defined by Eq. (11). Then as t— >iboo and 

x—>-^oo such that Aq ■.= \J —j> M and {x, t) G^ri ■={{x,t); x-4tA^ cos7„ :=/„(t) = 0(1)}, 



2 V t 

for those 7„e(f ,7r), 



where 
with 



Qix,t) = Qt,ix,t) + o{^^^), (13) 

g±(x,t) := Qg{x,t) + Qi{x,t) + Qf{x,t), (14) 

f^S ( f\ _ 2iAnsin(7n)cxp{^}exp{2i(A^(2tA^+;n(f)cos7n)+^^)} , . 

' ^ cosh(i^+2A2sin{7„)«n(i)-JJ) ' ^ ' 

^± = -iarga + arg5±(A„;Ao)+ E arg fcz^jlttSj) ' (16) 



2 

1&L± 



x± = -ln(A„sin7„) + ln|C„| -21n|,5±(A„;Ao)| +2^E l^ dtl^plti;!) ' ^^^^ 

5+(A;z) = exp{/,^ '^^^'^ - ^ ^T+i^f^ ^l. (18) 

5-{X;z) = exp{/r ^^f.^^}, (19) 

Qi{x,t) = y±^exp{i(,^±(Ao) + $±(Ao;t) + f )}, (20) 

z.(z) = -2^1n(l-|r(z)|2), (21) 

<P+{z) = i ln|£.2 _ z^\dln{l - \r{g)\^) - i J,^ ln\Q^ + z2|dln(l + \r{tg)\^), (22) 

r (^) = i /r H^-' - ^'|dln(l - \r{g)\'), (23) 
$±(Ao; t) = 4A^t T i^(Ao) ln|t| ± argr(ii/(Ao)) + argr(Ao) T 3i/(Ao) In 2 

+ (2±l)f + 2 E^arg (j^°:gj;:°aj) , (24) 



Qf (x, t) = - ^(^^^Ifl'^^l y±ig^{exp(zc^^(Ao; t)) + 2i cot(7n) cos(^^(Ao; t))}, (25) 

at ■= l9^|exp{iargg^}, 



bJl = |Cn||^^(An;Ao)|-^exp{-2A^sin(70;„(t)}exp 2E In t"!^ , (26) 



|A„-A,||A„+A,| 



±—r:,-rrrr< _9oT.rrA±('\ • \ „ _L 9 I ( An - A^) ( An + A£) 

(A„-AO(An+A;) 



argfifj = argC„ - 2arg5=^(An;Ao) + 2 E arg 



- 2A2(2iA2 +/„(i)cos7„), (27) 

^± ^ exp{ ^ } exp{2 A^ sin(7n )ln (t)-^ } , . 

2cosh(i:eL+2A2sm(7n);n(i)-JJ) ' ^ 

(/p±(Ao;t) :=arg5± + (/)±(Ao) + $^(Ao;t), (29) 
EzeL+ :=EzIn+i> EieL- ■=J2?=i > r(-) is t/ie gamma function |l|], and C(Ao) e5(M>Af;C), 
and, as i — > ±oo and x — > ±oo swc/i that fiQ := ^Jj > M and (x, i) G L5n := {(x, i); — x + 
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4tAlcosjn--=-ln{t) = 0{l)}, for those 7„e(0, f ), 

Q{x, t) = Q^{x, t) + o{^^^) , (30) 



where 
with 



X 



Q^ix,t) := Qi'{x,t) + Q''^{x,t) + Qf'{x,t), (31) 

nS'f^ t\ - 2A„sin(7„)exp{-%^}cxp{22(A2(2fA2+/„(t) cos7„)+0^O} , . 
^± ^' ~ sinh(^+2A2sin(7„)«„(t)+2S') ' 

0±' = -iargC„ + arg<5±(A;;/io)- E arg (g^^gg^) , (33) 

-ln(A„sin7„) + ln|C„| - 21n|5±(A;;Mo)| + 2^E l" ('|tlli!!t+^!) ' ^^^^ 

^+(A;.) = exp{/,^ .in(i+kMn ^ _ joo ^'ng-J^y) ^}, (35) 

<5,-(A;.) = exp{/r (36) 



gi'(a;, i) = -^^^ exp{i(<^±'(/.o) + $±'(mo; i) + tt)}, (37) 

Ki^) = -^ln(l + |r(iz)|2), (38) 

r'(^) = I /r - ^'|dln(l + HiQ)?), (40) 
$=^'(/xo;*) = =F J^(«Mo) ln|t| ± &xgV{iv{i^iQ)) + argr(i/io) =F 3i^(i/io) In 2 

-(2Tl)f-2 E-g(fe^5^ (41) 



Qf ^(x, = - 'I y^^|exp(iy^^(Mo; t)) -2i cot(7n) cos(y^^(/xo; t))}, (42) 

5n' := I5n'|exp{iarg5±'}, 

\9t'\ = |CJ|5±(A;;Mo)|-'exp{2A2 sin(7n)/n(t)}exp|2 E (43) 



|A„-A,||A„+Ai| 



(An-Az)(An + Ai) 



arg Qn = arg C„ - 2 arg 5P (A^; A^o) - 2 E arg 



2A2(2iA2 +/„(t)cos7„), (44) 



-±/ _ _ expi ^ 



exp{-%^} exp{-2A2 sin(7„)i„(t)-S±'} 



2 sinh( +2 A2 sin(7„ );„ (i ) 



(45) 



2 

(/.±'(/xo;t) := arg5±' + 4>^'{^lo) + $±'(/io;t), (46) 
and C(/xo)g<S(]R>m;C). 

Theorem 2.2 Let m(x,t;A) &e i/ie solution of the RH problem formulated in Lemma 2.1 
with the condition W\\c°=(^-^) < ^ '^'^^ Q{xji)> the solution of the DNLSE (Eq. (4)), be 
defined by Eq. (12) in terms of the function Q{x,t) given in Theorem 2.1. Then as >±oo 
and x^^oo such that Aq :=\^J —j> M and (x, t) G n„ := {(x, t); x— 4fA^cos7„:=Z„(t) = 
0(1)}, for those -in^{^,T^), 

q{x, t) = Qtix, t) exp{z arg g± (x, t)} + O (^i^al^) , (47) 
where Q^{x,t) are given in Theorem 2.1, Eqs. (14)-(29), 

arggi(a;,0 = -4 E 7i + 4arctan(77n|5n + cot7n) + 3^±(Ao) 

1&L± 
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+ 4, ± 



. '^(''^o) Ifln I sii(7n)(|gn P cos(.^^ (A() ;t) -7^ )-T;n sin(7n) cos((^^(Ao;t))) 



2Kt 



( (t/„ sin 7,1 + 1 (/J P cos 7„ ) 2 + 1 g J 1 4 sin^ 7„ ) 



2 j^c. VWm) (sft|i^±(o)} cos(e±(M; 0) + ^^{i?±(0)} sin(e±(/x; t))) (48) 



dz 



e 

d{r{z)\ 



z=0 



dz 



2=0/ ieL± J 



(50) 



e±(Ao;t) :=^HXo;t) + ^^{\o) + y±{\o), (51) 
with C(Ao) G »S(]R>m;'C), and, as t ^ ±00 and x — ±00 suc/t i/tai ^uq := 2\/f ^'^^ 
{x,t)eUn:={{x,ty, -x+4tAlcosjn---ln{t) = Oil)}, for those {0, ^) , 

q{x, t) = Qt:ix, t) exp{i arg q^{x, t)] + O (^I^ialfliim!) , 
where Q^{x,t) are given in Theorem 2.1, Eqs. (31)-(46), 

COt7n) +y'±il^o) 



(52) 



argg±'(x,t) = 4 E 7i + 4arctan(?7„|5^ 



,±/l-2 



ieL± 



i^Cw) Ifln'l s™(7n)(|gn cos(ip^' (^Q ;f )+7n )+)?n sm(7n) cos(y^' (/JQ 



„±'|2 



((r/„ sin7n-|gJ'P cos 7„)2+|gJ'|4 sm27„) 



±1 C "^^^ (^{^^'(0)} cos(e±'(^; 0) + 9{ii±'(0)} sm(e±'(^; t))) ^, (53) 



i2±'(0):= 



dz 



2=0 



d2 



lexp^ 

2=0/ 



-4i E 7« 



G±'(mo;0 := $^'(/io;i) + </'^'(/^o) +:^i(/^o), 



(55) 
(56) 



and C(/xo)g5(E>m;C). 



Theorem 2.3 Let m{x,t;X) be the solution of the RH problem formulated in Lemma 2.1 
with the condition lkll£oo(R.C) < 1 '^f^d u{x,t), the solution of the MNLSE (Eq. (2)), be 
defined by Eq. (10) in terms of the function q{x,t) given in Theorem 2.2. Then as i— ^±00 



and X— *-±oo such that := y ^(f — 7) >M, f > ^, sgM>o, and (x, t) Gr2„ ■={ix, t); —x + 
t(i-2A2 cos jn) ■.=Tn{t) = 0{l)}, for those jn G (f , tt), 

U{X, t) = Vtix, t)wt (X, t)+0 (^ 0(Xo)(}n\t\)- ^ ^ ^57^ 



where 
with 



v^{x,t) ■.= vi{x,t)+vg{x,t)+vg^{x,t), (58) 

yS( 4.\ ^ V2iA„sin(7„)exp{^}exp{2i(A2(tA2_^r„(f)cos7n)+0^)} , . 

^± = - i arg Cn + arg <5± ( A, ; Ao) + ^ E^arg ([tl^jit+I;)) ' (60) 



± = -ln(A„sin7„) +ln|C„| -21n|5±(A„;Ao)| + 2 E In ( ) ' (61) 



|A„-A;||A„+A,| 



vg{x,t) = y±^exp{i(</.±(Ao) + $±(Ao;t) + f )}, (62) 
$±(Ao; t) = 2X^t T i'(Ao) ln|t| ± arg r(zz/(Ao)) + arg r(Ao) T 2i^(Ao) In 2 
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+ (2±l)f + 2 Earg(ga:iM^), (63) 



^c^^^i) ^ _li^J_Md^±^|exp(i^±(Ao;t)) + 2zcot(7n)cos(^±(Ao;t))}, (64) 

5^ := Iffn |exp{iarg^^}, 

\gt\ = \Cn\\S^{XnrXor'eM-'^^lMln)Tnmewl^E^^ (65) 



arg gt = arg C„ - 2 arg 6^ (A^; Aq) + 2 arg 

leL± 



(An-Aj_)(An+A£) 

(A„-AO(An+Ai) 



- 2A2(tA2 +^„(t)cos7„), (66) 

^± _ exp{^}exp{2A^ sin(7n)Tn(t)-a^^} /g^N 

ipt(Xo;t) :=arg5± + ,/.±(Ao) + $±(Ao;t), (68) 
u;±(x,t) = expii -4 7z + 4arctan(r?„|3±|-2 + cot 7„) + 3^±(Ao) + 2jj(4A§s 



_!_]_) _|_4 /_[_ KAq) |g^| sin(7n)(|gn P cos(yj^(Ao;t)-7n)-T;^sin(7„) cos(y?^ (Ao;t))) 
Y Agt ((»?„ sin7„+|5*|2cos7„)2+|gJ|4sm27„) 

e±(Ao;t) :=$±(Ao;t) + </.±(Ao)+3^±(Ao), (70) 

and C(Ao) G 5(IR>m;C), and, as t^itcx) and x^^oo or ±oo such that ^^0'-= \l \{\ — j)> 

M, f <i, sgM>o, and (x, t) G 0^ := {(a;, t); x - t(i - 2A2 cos 7„) :=-!„(*) = 0(1)}, /or 
those 7nG(0, |), 

«(X, i) = Vtiix, t)wti{x, t) + 0( ^(^o)(ln|^|)2 -| ^ (71) 

^;i'(x, t) := 4'(x, t) + 4'(x, t) + vf'{x, t), (72) 

+\ _ V2A„ sin(7„) exp{- cxp{2i(A;^^ {tAl+?„{t) cos 7n)+0^O} /^oN 



where 
with 



yisinh(i^+2A2 sin(7„)Z„(t)+2^J') 

= -i arg C7„ + arg S^(K; fio) - ^£arg(g^ESfe$S) , (74) 
xf = -HAnSmjn) + ln\Cn\-2\n\d^{K;no)\ + 2Eln(^^^^ (75) 

$=^'(/io;*) = 2p,Qt^i'{ip,o)ln\t\ ±argr(izy(i/io)) +argr(i/io)=F2^(i/io)ln2 

- (2 T l)f - 2 E argffeSISj) , (77) 



4 ,g^^^V(w-AO(M+AO 



^5C.(^^^)^_MH;0M:^^^ilgl{exp(i^±^ (78) 

9n ■= I5n'|exp{zarg5±'}, 



l^ri = |gn||.5f(A„;/2o)|-^exp{2A^sin(7n)/n(t)}exp^2 Eln( j;;-gjj^;;;gj l^ (79) 
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(An-Aj)(A„+Ai) 



arg gt:' = arg C„ - 2 arg (A^; /io) - 2 E arg 



- 2A2(tA2 +Z„(t)cos7„), (80) 

_ exp{-i^} exp{-2A2 sin(7„)?„(f)-x±'} 
2sinh(i^+2A2 sin(7„)«„(t)+J±') ' 

^±'(//o;t) :=argff±' + ,/)±'(/io) + $^'(/io;t), (82) 
w;i'(x,i) = expii 4 E 7i + 4arctan(r7„|^^'|-^ - cot7„) + 3^^t(/io) + 2|7(-4/xgs 



-1-11—4 /-[- '"(jf'o) sm(7n)(|3n cos(y^'(w;*)+7n)+??n siii(7„) cos(y^'(/Jo;t))) 
V^^^* ((»?nSin7„-|flrpcos7„)2+|s^rsin27„) 

(3^{^^'(0)} cos(e±'(/.; t)) + Q{i?±'(0)} sm(e±'(/x; i))) 1^)} , (83) 

e±'(/2o;t) := ^^'(/Io;i) + </'^'(mo) + 3^i(/io), (84) 
and C(/2o)g5(]R>m;C). 

One possible application of the asymptotic results obtained in Theorems 2.1-2.3 is 
associated with the so-called "soliton scattering", namely, the calculation of the position 
and phase shifts of the rath soliton (1 < n < A^) for Q{x, t), q{x, t), and u{x, t) in the presence 
of the continuous spectrum: other physical applications of these asymptotic results include, 
for example, the calculation of the temporal and spectral intensities for the solutions of the 
DNLSE and MNLSE. 

Corollary 2.1 

(A) Q{x,t): 

Axf := i2nn)-Hx+-x-) 

N 



T sgnd - n) In - Xi\\Xn + Xi\ \ ^XuM 



A</.f := 2(0+ 

N 



2 Vssna nUrJ^^^i^^M^^UarJ^^^^ 
2 <; ^ sgn(^ - n) argi ^ _ ^^^^^^ ^ ^^1 +argU_^^^.^^, 



1=1 



Ax^ := {2r^n)-\xt' - X-') 

= v-n' I E sgn(i - n) mfl^- - ^11^- + 'J] -mf ^"^'^ 

A0f' := 2(^+'-^-') 



-2^ 



1=1 



>; 



(B) q{x,t) (DNLSE): 



Axl =Ax^ , 
A</)^ =A<f>^ -4j:sgn{l-nhi+y+{Xo)-y-{Xo), 



1=1 
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1=1 

(C) (MNLSE): 



Axf 


= Axf 1 


Ao- 


-Ao' 


A0f 


= A<^f 1 


Ao- 


»Ao' 


Axf 


= A.t9^'| 












AC" 


= A^f 







Proof. Follows from the definition of soliton position and phase shifts given in Q and 
Theorems 2.1-2.3, Eqs. (16), (17), (33), (34), (48), (53), (60), (61), (69), (74), (75) and 
(83). ■ 

Remark 2.5 The expressions for the soliton phase shifts given in Corollary 2.1, namely, 
A(/)^^, A(/)^'^', A(f)f , A(j)i^\ ^(l>^ ^ and A(j)'f\ 1 <n<iV, are to be understood mod(27r). 

Remark 2.6 For the asymptotics of the C-valued functions Q{x,t), q{x,t), and u{x,t), 
one must actually consider four different cases, depending, respectively, on the quadrant of 
the (x, t)-plane. In this paper, the proof of the asymptotic expansions for Q{x, t) and q{x, t) 
(resp. u{x,t)) is presented for the cases (x, f) — > (=Foo, ±oo) (resp. (x, t) ^ (±oo, ±oo)) such 
that Ao > M and (x, t) G 0„ (resp. \o> M and (x, t) S Cln) for those 7„ S (f ,7r): the results for 
the remaining domains of the (x,t)-plane are obtained analogously. If the conditions on 7„ 
stated in Theorems 2.1-2.3 are violated, then (x, t) G {]R^\il„, ]R^\J^„, M^\On, but 
the asymptotic results stated still remain valid although the second terms on the right-hand 
sides of the asymptotic expansions become the leading-order terms of the corresponding 
asymptotic expansions, while the remaining terms are exponentially small and negligible 
with respect to the given error estimations. 



3 The Model RH Problem 

In order to simplify the asymptotic analysis of the original RH problem formulated in 
Lemma 2.1, a simpler, model RH problem (see Lemma 3.3) is derived in this section. As an 
intermediate step towards the formulation of the model RH problem, it will be convenient 
to derive an "extended" RH problem (see Lemma 3.2): the general idea pertaining to 
the transformations from the original RH problem to the model one is elucidated in the 
paragraph following Lemma 3.1 (see below). 

Proposition 3.1 (Q) In the solitonless sector [Z^ = 0), as t ^ +oo and x — oo such 
that Ao:=^y^>M, 

m(x,t;A) = A(A) + 0(^), 

where A(A) := (5+(A;Ao))''^ 

5+(A;Ao) = ((^)(^))"exp{ E (p^A) + P/(A))}, 

P±{^)-2riJ0 ^^[l~\r{Xo)[^ )T^' P±K^) - J±ioo (?rA) 2H' 

u := u{Xo) is given by Eq. (21), ll('^^(-;Ao))'^^||^oo(C;C) sup^g(j-< |((5+(A;Ao))=^^| < oo, 

^5+(ibA;Ao)^ =(5^(A;Ao), the principal branch of the logarithmic function is taken, ln(/i— 
A):=ln|/i-A|+iarg(/i-A), arg(/i- A) e (-vr, vr), anc? C(Ao) e5(M>M;Af2(C)). 
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Remark 3.1 For notational convenience, until the end of Sec. 5, all explicit x,t depen- 
dences are suppressed, except where absolutely necessary, and (5+(A;Ao) :=(5(A). 



Lemma 3.1 There exists a unique solution m^(A): C\(2^rfUr) ^SL(2,C) of the following 
RH problem, 

1. m^(A) is meromorphic y X(zC\T, 

2. 

m^{X) = m^{X)v^{X), A G f , 

where 



^A^^ g-^e(A)ad(.3)[ r{X)5^iX)6^{X) \ 

\ -r(A)(5_(A)r^(5+(A))-i i6^iX))-H^iX) J' 

and 0{X):=X'^x + 2XH, 

3. m^(A) has simple poles at {itAj, ±A.j}^;^ with {l<i<N) 

res(m^(A);Ai) = lim m^(A)wi(5(Ai))~V_, 

res(m^(A);-Ai) = -cT3res(m'^(A);Ai)cj3, 
res(m^(A);A;) = linLm^(A)W(<J(\))V+, 

res(m^(A);-Aj ) = -(T3res(m^(A);Aj )(T3, 

4. as A ^oo, AG C\(ZrfUf), 

m^(A) =I + 0(A-i); 

moreover, Q{x,t)=2i lim (Am^(x, t;A))i2 is equal to Q{x,t) in Lemma 2.2, Eq. (11). 

A— »oo 

Proof. Let m[X) be the solution of the RH problem formulated in Lemma 2.1. Define 
m^(A): = m(A)(A(A))-i. ■ 

In order to motivate Proposition 3.2 and Lemma 3.2 (see below), consider the trajectory 
of the nth soliton with 7„ G (-1,71") in the (x, t)-plane which belongs to the set Vt^ := {(x, t); x— 
4tA^ cos 7„ = 0(1)}, and note from Lemma 2.1 and the soliton ordering in Sec. 2 that, as 
t — > +00 and X —CO such that Xq> M and (x, t) G ^n'- (1) '^{vi\n„) ~ C(exp{— 8t7/j(^j — 
C„)})^OVi<n (iG{l,2,...,n-l}); (2) 3f?(t;,|oJ ^ 00 Vi > n (i G {n+ 1, n + 2, . . . , iV}); 
and (3) 3f?(vj|r2„) ~ for i = n. Thus, along the trajectory of the arbitrarily fixed nth 
soliton, there are exponentially growing polar conditions for solitons i with n + l < i < A^. 
One must effectively deal with such growing polar conditions in a self-consistent manner. 
In a recent paper [pO| devoted to the asymptotics of the Toda rarefaction problem, Deift et 
al. showed how this could be done: they noticed that it is possible to replace the poles with 
the exponentially growing polar conditions by jump matrices on small, mutually disjoint 
(and disjoint with respect to T) circles such that these jump matrices behave like I -|- 
exponentially decreasing terms as t^+00. Thus, instead of the original RH problem, one 
gets a new, "extended" RH problem with 4(A^ — n) fewer poles, and 4(A^ — n) additional 
circles with jump conditions stated on them. Finally, by removing the added circles from 
the specification of the extended RH problem, one arrives at the model RH problem: the 
estimation of the "difference" between the extended and model RH problems shows that 
the solution of the model RH problem approximates the solution of the original one modulo 
terms which are decaying exponentially as t^+oo. 

Proposition 3.2 Introduce arbitrarily small, clockwise- and counter- clockwise- oriented, 
mutually disjoint (and disjoint with respect to T) circles Kf and Lf, n-\-l<j <N , around 
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the eigenvalues {±Xj}jL^_^_i and {±\j}jL^^i, respectively, and define 



N 



m^(A), AGC\(ru( U {KfuLf))), 



m 



^(A)(l- ^ 



i=n+l 



(AtA.i A€intKf, n + l<!<W, (85) 



AeintLf, n + l<i<N. 



Thenm^iX) solves a RH problem on (cri;\U,^„_^i({±Ai}U{±Ai}))U(U,^„_^i(iv:^ULf )) with 
the same jumps as m^(A) on T, m^_(A) = m^(A)v^(A), and 



m^(A) 




(A)(l+^ 



(A=FAi) 



XeKf, n + l<i<N, 
XeLf, n + l<i<N. 



Proof. Follows from Lemma 3.1 and the definition of m''(A). ■ 

Remark 3.2 The superscripts it on {Kf}fL^_^i and {Lf}fLn+ii which are related with 
{±Ai}^„^i and {±Ai}^„_,_^, respectively, should not be confused with the subscripts ± 

appearing in the various RH problems in Sees. 3-5, namely, m±(A), m^(A), mj_(A), m?|_(A), 
X±(A), ^±(A), and xl(A). 



Remark 3.3 Even though the exponentially growing polar (residue) conditions have been 
replaced by jump matrices, it should be noted that, along the trajectory of soliton n, these 
jump matrices are also exponentially growing as t — > +oo. These lower/upper diagonal, 
exponentially growing jump matrices are now replaced, through a sequence oi N—n sim- 
ilar transformations, by upper/lower diagonal jump matrices which converge, along the 
trajectory of soliton n, to I as t^+oo. 



Lemma 3.2 Set 



mtt(A):=< 



N 



N 



m\X) n (d;+(A))-'^3, AGC\(ru( U (i^fuLf))), 

l=n+l i=n+\ 

m\X){Jj^±{X))-^ n {duWy^, AGintKf, n + l<i<N, (86) 



l=n+l 
N 



m\X){Jj^±{X))-^ n (^^^_(A))-'"^ A G intL^ , n + 1 < i < iV, 



l=n+l 



where 



N 



di^{X):=^^^^^,, AeC\( U {KfuLf)l n + l<l<N, 



i=n+l 
N 



^^"(AlAy'^ AG U inti^^ n + l<l<N, 



i=n+\ 
N 



(87) 



^^^^'^ , , AG U intLf, n + l<l<N, 



=n+l 

and the SL{2,<C)-valued, holomorphic in inti^^^ and intL^, respectively, functions J^±(A) 

i 

and J^±{X), n+l<i<N , are given by 



J^±(A): 



l = 7l-\-l 



(A=FAi) 
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' (A^A,) n ^ 

i=n+l , i'^i 



(CT))^ dr;w 



with 



w(5(±A,))2 n ^ 

11 (d, (A))-^ 11 7^ 
l=n+l ^ — 1=71+1 


(A=FAi) 



N 



c} = {vi)-HS{±Xi)f{di_{±Xi)f n (^^/+(±A,))^ n + l<z<iV. 

l=n+l 



(88) 



T/ienmtt(A):C\((^rf\uiI„+i({±A,}U{±A,}))U(fu(uiI„+i(i^fuLf))))^SL(2,C) soZ^es 
the following, extended RH problem on {a£\uf^^^^{{±Xi}U{±Xi}))U{ufL^^;^{KfuLf)), 

w;/tere 



t;«(A)|p: 



(l-KA)r(A))|^ 



r(A) 



iV 

n 



\^ d- (A)d+ (A) ^^1^1_^1^(A- A, )( A+ Ai ) 



(A-Ai)(A+Ai)^ ^ 



•5+ (A) 

5^Ay 





-A,) 


(A+Ai)\ 




-Ai) 


(A+AOi 



.«(A): 



1 + 



1 + 



(a=fa<) 



■2\2 N 



2\i 



n 



aT'-a? 



I=n+1 



N 



(A=FAi) 

with polar (residue) conditions, 



(v-i)-H5i±_x,))-^ fxj-x 

2Xi 



a+, AG U 

i=n+l 



n 

=n+ 



!=n+l VA; —A, 



N 



(89) 



AG U ^ 



± 



i=n+l 



N / \ 2 

res{rni{X);X,)=l\mm\X)vMX^))-' ^n^^i^^^ a-, 1 < i < n, 



(Ai-Ai)(Ai+A;) 1 l<i<n 
-o-3res(m''(A);Aj )as, 1 < ^ < n, 



res(m''(A);— Aj) = — (T3res(m''(A);Aj)(T3, 1 < i < n, 

_ , _ Af /- - - ,-\ 2 

res(m«(A);AO= linL"^*(A)W('5(A0)' H 

A— +Ai i=n+l ' 

res(m''(A);— Aj ) 

and, asA^oo, AGC\((^d\uiI„+i({±A.}U{±Al})) U (ru(uiI„+i(i^±UL±)))), 

m«(A) = I + 0(A-i); 

moreover, Q{x,t) = 2i lim (Am^(x, i;A))i2 is equal to Q{x,t) in Lemma 2.2, Eq. (11). 

A— >oo 

Proof. The proof is presented for the eigenvalues {Aj}^„_|_^, around which arc defined the 
small, clockwise-oriented, mutually disjoint circles the proof for the eigenvalues 

{~Aj}^„_l_^ and {±Ai}^„_|_]^ follows in an analogous manner. From the definition of m^{X) 
and Proposition 3.2, one sees that, on {-^^j^liln+i' "^''(A) solves the following RH problem 
(AGUiI„+ii^+), 



N 



m: 



■+ 



(A)=m«_(A) n {di_{X)rJ^^{X)[l + 



l=n+l 



V^{5{Xi))- 

(A - Xi) 



N 



-a\ n {di^iX))- 



-<T3 



l=n+l 



jump matrix 
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Demanding that the above "jump matrix" be equal to the following upper triangular form, 
I+p^z^i^+j n+l<i<N, one shows that 



' vm^i))-'cj n ^ n Ttv^ ^ 

c'i-W 'ii-W {x-Xiy^di_{x)di^{x) {x-Xi) di_{x)di^{x) 

(^-^*) ^ (=y+i ^ ^ .=y+i ^ 



VT^ i=n-|--L T ■ " -r ■ ■ i=n-\-i ^ 

Note that det(J^+(A)) = 1 (n + 1 < i < N). Defining, for n + 1 < Z < N, di JX) and di_{X) 

i 

as in Eqs. (87), and choosing cf, n + 1 < i < N, as in Eq. (88) (with one gets the 

expression for J^+(A) (which is holomorphic VA G UfLj^.iintK^) given in the Lemma; 

i 

also, because of the symmetry properties of (5(A) (Proposition 3.1), Cf = {vi)~^{S{±Xi))~^ 

N 

• {di_ (zbAj))"^ n {di+ (=t Aj))~^. The remainder of the proof is a consequence of Lemma 3.1, 

Proposition 3.2, and the definition of m\X). ■ 

Remark 3.4 Even though, along the trajectory of soliton n, all the initial, exponentially 
growing nilpotent residue matrices have been replaced by jump matrices which tend to I 
as t^+oo, i.e., 3e gR>o such that VzG {n+1, n+2, . . . , A''} |(vi|n„)~^| ~ C){exp{—et}), it 
does not necessarily follow that elements in the solution of the extended RH problem for 
m^{X) cannot grow exponentially; for example, note that the (2 l)-elements of J^±(A) and 

i 

the (1 2)-elements of Jj^±{X), n+l<i<N, grow exponentially. 

i 

By estimating the error, along the trajectory of soliton n (1 < n < N), when the jump 
matrices on {K- ,Lf}fL^_^_i are removed from the specification of the RH problem for 
m'^{X), one gets the following — asymptotically solvable — model RH problem. 

Lemma 3.3 Let x(A) solve the following RH problem on cJi;\U^„_,_i({±Aj}U{±Ai}), 

X+(A) = x-(A)e-*^W-'i('^3)^tt(A)|p, A G f, 
with polar ( residue ) conditions, 

res(x(A);-Ai) = -(73res(x(A);Ai)(j3, 1 < i < n, 

j\f / \2 

res(x(A);Al)=liin_x(A)W((5(Al))2_n (^^^^^^\ a+, l<i<n, 

res(x(A);-Ai ) = -cr3res(x(A);Ai)cr3, 1 < i < n, 
and, as X^oo, AGC\(Uf=i({±Ai}U{±Al})uf ), 

X(A) = I + 0(A-1). 

Then as t ^ +00 and x —00 such that Xq> M and {x,t) € the function E[X) := 
m''(A)(x(A))~-^ has the following asymptotics, 

E{X) = I + C'(F(A;Ao) cxp{-a6t}), (90) 

where \\F{-,Xo)\\^^^c.j^^^q^ < 00, 11^(^5 •)II/;-(M>m;M2(C)) < ^(A;Ao)~o(^) as 
A ^ 00 with C(Ao) G C°^{R>m;M2{C)), a := 8min{r?i}^^„_^i (> 0), and b := min{|^„- 

^i|}i^n+l- 
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Proof. Writing, for n + l<i<N , Eqs. (89) in the following form, 

N 

I+(W)-^W^±(A)a_, AG U Lf, 

» i=n+l 

consider the "error function" E[\) defined in the Lemma. One notes that: (1) det(S(A)) = 1; 
(2) E[\) has no poles; and (3) E{\) solves the following RH problem on the oriented contour 
Ss:=uiI„+i(i^+Ui^ruL+ULr), 

(W)-'W,.(A)(-£)-W )), A . Lf. 

n+l<i<N, and, as A^oo, X£C\T.e, E{X) = 1+0{X~^). Now, writing the RH problem 
for E{X) on the oriented contour S^; in terms of an equivalent system of linear singular 
integral equations, using the explicit asymptotic solution of the model RH problem for x(A) 
given in Sec. 4, recalling that, as t — > +oo and x — > — oo such that Aq > M and (x, t) G ^n, 
(^i|n„)~^ ~ C(exp{— 8t7/i|^ri— ^^+1 < ^ < -^j and proceeding as in the proof of Lemma 3.3 
in IQ, one deduces the estimate in Eq. (90). ■ 



4 Asymptotic Solution of the Model RH Problem 

In this section, the asymptotic (as t — > +oo and x/t ~ 0{1)) solution of the model RH 
problem (Lemma 3.3) for the Schwartz class of non-reflectionless generic potentials (r(A) G 
5(r;C)) is presented. Before doing so, however, recall the following well-known fact from 
matrix RH theory |Tl| , 

Proposition 4.1 The solution of the model RH problem (Lemma 3.3), x(A):C \ (L U 
(U"^]^({±Aj}U{±Aj}))) ^SL(2,C), has the following representation, 

, /■ x-(g)(^»(g)lf-i) de 



r 



where 



, , M \ _ T _|_ ( res(x(A);Ai) g-3res(x(A);Ai)q-3 res(x(A);Ai) g-3res(x(A);A,)o-3 \ /Qr,\ 

Xd[X)-i + ^^y-j^ ^+A7^ — + ^3^17^ — )' ^ ' 

The solution of Eq. (91) can he written as the following ordered product, 

x(A) = Xd(A)x^(A), (93) 
where Xd{X) is given by Eq. (92), and X^{X) solves the following RH problem: (1) x'^{X) is 
piecewise holomorphic VAGC\r; (2) X+(A) = X-(A) exp{— i^(A)ad((T3)}(u^(A)|p), A G F; 
and (3) as X^oo, AGC\f, x^(A) = I+0(A-i). 

Remark 4.1 From Proposition 4.1, Eq. (93), it is seen that, in order to solve the model RH 
problem, explicit knowledge of Xrf(A) and x^(A) is necessary. The determination of x^(A) is 
technically the more complicated of the two: actually, the determination of Xrf(A) depends 
on the explicit knowledge of x'^(A) (see Proposition 4.2); hence, the asymptotic solution of 
X'^(A) is presented first (see Lemma 4.1). 

In order to more fully comprehend certain elements of the proof of Lemma 4.1 given 
below, the Beals-Coifman |l^ formulation for the solution of a (matrix) RH problem on 
an oriented contour is requisite: a self-contained synopsis of this formulation as it ap- 
plies to the solution of the RH problem for x'^(A) stated in Proposition 4.1 now follows. 
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Writing the jump matrix in the following factorized form, ?;'*(A)|j^ := (I — tf^^^(A)) ^(1 + 
^XtW), Aef, where ti;^^t(A) G nfcg{2,oo}^^(r;M2(C)) (with \\wt,ti-)\\r,^^^^ ^^cHf-,M2{C))-= 
Efce{2,oo} lkx,t(-)ll£fc(f;M2(C)))' respectively, are nilpotent off-diagonal upper/lower trian- 
gular matrices, define Wx,tW := w~ ^{X) + ^{X) , and introduce the operator C^^^ on 
£2(f;M2(C)) as C^^J:=C+ifw-t) + CMwXt), where f G C\f;M2iC)), and C±: C\f; 
M2(C))^£2(p.^2(C)) denote the Cauchy operators, (C±/)(A) := lim (^^• 



A' e± side of? 



Theorem 4.1 ([0) ///u'^(A) Gl©£^(r;M2(C)) solves the following linear singular integral 
equation, 

(Id-C^.,K = I, 

where Id is i/ie identity operator on I©£^(r;M2(C)), t/ien i/ie solution of the RH problem 
forx^W is 

x(A) = i+/ — TT — — , AeL\r, 

where ,i%X) = x%{X){l+wltWr'=xUm-w~M)r^- 

Lemma 4.1 Let eo denote an arbitrarily fixed, sufficiently small positive real number. For 
G {0, ±Ao}, set AA(H;eo) := {A; |A — H| < eo}. Then as t ^ +oo and x — > — oo such that 
Ao>M anc? A G C\Ui-<g|o,±Ao}-^(^;fo); X'^W has the following asymptotic expansion. 



= I + i/^(x4l^ + X^)(-P{-(^-(Ao) + $-(Ao;t))K 
+ exp{i(0+(Ao) + $+(Ao;i))W) +0^^^^'^°)^"*^ 



t 

where i^(Ao), (/>"'"(Ao), and <^"'"(Ao;t) are given in Theorem 2.1, Eqs. (21), (22), and (24), 
l|G('^o)||^.(C\u.,,o,.Ao>^(H.o);M.(C))<°«' G(A;.)e5(M>M;M2(C)), G(A;Ao)~o(^) 
as A ^ oo with C(Ao) G 5(M>m;-^2(C)), and satisfies the following involutions, x'^{~^) = 
cr3X''(A)(T3 and x''(A) = f7ix'=(A)cri. 

Proof. In Sees. 5 and 6 of 0, it was shown that, for A G C \ Ui-(g|o,±A(j}-A^(^;eo)i as 
t^+oo and x^— oo such that Ao>M, for arbitrary /'gZ>i, 

mseMs ai (ma) 

where: (1) G M^. denotes the set of vectors with k components each of which take the 
values ±1 (card(Mfc) = 2^^); (2) ai{rhz) = miAo -|-m2eexp{ — ^^^}, au{rh^) = miXo, and 
b{fh2) = migexpj ^^^^ }, where e is an arbitrarily fixed, sufficiently small positive real 
number; (3) := /U^(?)Il'> — /"'(?)ll7> with L' = {A; A = 2exp{f }, u G 

(— e, e)} U (U;g|-|-i}.{A; A = IXq + uexp{ — ^}, u G (— and where /x'^(-) is the solu- 
tion of the Beals-Coifman [|l3| linear singular integral equation (Theorem 4.1); (4) 0(?) = 
2,2(,2_2A2); (5) 5(,) = ((Mii)(i^))-exp{E,e{±}(pK?) + ft(0)}, ^ := KAo), /^±(?) = 

1 r=t^01nridlMl^-^ and n, ^r^ - ln(l-r(g)r(g)) , /„n t-, -1,-1,1 / x _ 1,1,1 / x _ 

2¥iJo in(^jq^:p^J^^, and pil^g - J^i^ ^^^^^ K'^) \S) - K- \S) - 

-r{s)V{s), n-^'''\,)=n^^~^'\,) = 7^-l'-l'-l(<^) =^^'^'-H<?) = (T^i'^'H?))*, 
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and 7^-l.l'-l(0 = 7^l'-l.-l(<^) = (7^l'-l'l(?))^ where V{z) := nil,+i ({f^^^jg^) and 
a(-) = means that a(-) is the same piecewise-rational function as except with 

the complex conjugated coefficients; (7) TZ^'^{<^) = 7^"^'^(?) = ([l'^^.'^^^^) - (i+'fr^^^p) and 
7^l'-l(^)=7^-l.-l(^) = -(7^l'l((^))*; and (8) c(Ao)G/:°°(M>m;M2(C)). Since, in the above 
expression for x'^W: the estimation of all the integrals is analogous, without loss of gener- 
ality, the following integral is considered, 

Jo := / A(?)^oWii, 

_ ITT 

where :=/""(OIl', and Bo{<i) := - ^'^'^ «^p{-2»f%)}Kg)y W Begin by estimating Bo{^) 

"im ' expand So(?) m a Taylor series about Aq; (2) make the following change of variable 
0, q:=q{w) = XQ+w{lQ\lt)-^/'^ , and express the expansion obtained in (1) above in terms 
of w; and (3) use the following identity, ah={a—l){h—l) + {a—l) + {b—l)+ 1. Carrying out 
steps (l)-(3), one gets that. 



r4eL4 ^==1 



(A-Ao) (A-Ao)2J ^3/2+^„ 



where £ denotes the set of vectors with k components each of which take the val- 
ues and 1 (card(Lfc) = 2^=), s{Xo) := Eie{±}(ft(Ao) +P/(Ao)), 7^(Ao) := -r(Ao)P(Ao), 

n'iXo) = -r'(Ao)P(Ao)-r(Ao)P(Ao) Etn^i{ '''X"-tj''^^ + ^tSSt^}' ^ 



|z=Ao) 



:= exp{2Afc(Ao + ^7=p)} " 1> {1, 2}, 



'4 



and Cl'(Ao) G 5(M>m;C), i S {1,2}. Now, proceed to estimate /i'^(?)|L' for <; G (Ao,Ao + 
eexp{— ^}): for this, the Beals-Coifman [l^ formulation for the solution of a RH problem 
on an oriented contour is necessary (Theorem 4.1 and the paragraph preceding it for discus- 
sion and notation); in particular, one has to estimate the functions (<j) on L'uL'. In Sec. 5 
of |], it was shown that, on (Aq, Ao+eexp{-f }), = -(<5(?))^'^('^3) exp{-it6'(<;)ad(o-3)} 

• r(<j)P(<j)(T+ and =0; hence, from the Beals-Coifman ||l^ formulation, for any /G 

£^(L';M2(C)), Cut^tf = C^{f w^f-). To estimate w^^{(^), one proceeds as follows: (1) recall- 
ing that r(0G5(f;C) and ||r||^^(K.Q < 1, expand iS{<i))^'^^''^\ for ? e (Aq, Ao+e exp{-f }), 

via an integration by parts argument; (2) expand exp{—it6{q)ad{a3)}r{<^)T'{<^) in a Taylor 
series about Aq; and (3) change variables [^, <^ := <^{w). Carrying out steps (l)-(3), one 
shows that. 
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where (/<(Ao) = 21. InAo + ^/o"" ln\z'-Xl\dH'i^-\r{z)f)-^J^ln\z' + Xl\dln{l + \r{iz)\^), 
vqo^WjXq) :=7^(Ao), vio{w;Xo),vii{w;Xo), and V22iw;Xo) are nilpotent matrix polynomials 
whose elements are sums of products of terms of the type and (Inw)^, j gZ>i, k^Ti^Q, 
with Ao-dependent coefficients, and, for 0<j<2, 0<k<j, 

1 1 exp{-i(-)2} exp{2izy ln{-)}vjki-,Xo)\\^^^^^ ^ ^j£i(L'A{o};C) < 

"^^^^ "^■)"n,e{i,2,^}/:Hi^\{o};C) •= ^ie{i,2,oo}\\i-)\\ciiL'M0hCy ^^^^^^ 
scaled and shifted version of -L'\{Ao}. Hence, for any /g£'(L'^;M2(C)), /g{2,cxd}, 

(C.,,J)(^5;Ao) = E E (^(Cy)(u;;Ao) + o(^(C^2/)(t^;Ao)), 

^^^'^^ exp{2i</.(Ao)}'^Sfc.''''^^''^0'' Jl' {z-X') 27ri ' ^ - - ^' ^ - 

A' side of lJ. 

A;<j, and l|Cj-fc(-;Ao)||;vi(L',\{o};A/2(C)) -''^i('^o) '^^^'^ M{»;M2{C)) denoting the space 
of bounded linear operators acting from £'(»;M2(C)) into £^(»;M2(C)), /g{2,oo}. Accord- 
ing to Theorem 4.1, ^'^(•) satisfies the following linear singular integral equation on L'UL', 
( Id — Cuj^j.)fi'^ = I; hence, ^'^ = (Id — Clu^ ^)~H. It was shown in that, as t ^ +oo and 

x^-oo such that Ao>M, ker(Id - C'^„,,t) = and ||(Id - C'5o(-;Ao))~"^||;vj{L(,\{o};A/2{C)) - 
^2(Ao) < oo. Using the method of successive approximations, one shows that, as t^+oo 
and X— > — oo such that Aq > M, ^'^{■) can be expanded in the following Neumann-type series 
(see, also, Part II of |l| , and |l| ) , 



where ^go(^^;Ao) := (Id - cIq{wM)) ^I, ||(Id - C^o(-;Ao)) I£2(l^\{o};A/2(C)) < 
ll/"ifc(-;Ao)||£2(i/\|o};Af2(C)) < oo, 1 < J < 2, < < j: an explicit expression for (Id - 
Cqq{'w;Xq))~'^1 in terms of parabolic-cylinder functions was given in Sec. 7 of (see below). 
Making one more change of variable, q = \/2'5exp{^}, and recalling the definition of /q, 
one shows that 



where 



-^0 — h/2 — Io,a + Io,b + Io,c + Io,d + £r 



h/2 ■■= Ya{X,Xo;t) fiU^MQ^'^e-^'^Wdg, 



f ^ _ ITT 4 _ 

Io,a:=YaiX,Xo;t) E Iof^mi'^MUiPki^)y'0'''e-^y^a+de, 

^ _ ilV 4 _ ITT 

Io,,:=n(A,Ao;t) E /oV8o(^;Ao) n(Pfc(^)y^^>'^^+^e-^'/Vde, 

r46L4 ^ ^=1 

^ _ ITT _ ZTT 

/o,c := yc(A,Ao; t) E r{/^fo(^;Ao) + /ifi(^;Ao) Int} 

/4G]Lf4 

4 _ ITT 

X n(PA.(^))''=^'2*"e-^''/V+df?, 
k=l 

^ _ ITT 4 _ ilV 



+ o(!pS{E^^^^}/o"{Mfo(^;Ao) +/^fi(^;Ao)lnt},2-+ie-^V2^^rf^ 
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S := (32e^Aot)^''^, the prime on the summation in the expression for /o,a means that the 
term corresponding to (/i, ^3; ^4) = (0, 0, 0, 0) is omitted from the sum, 

Ya{A,Ao,t)- (^_^^^)^32^., , Yb[A,Ao,t)- (A-Ao)(32A^t) ' 

V(\ X^-f-^l — ?^(^0'*)g~^^(^o) v.a X •f-'i — t^;(Ao;f)7e(Ao) 
^^^-^'^0' ^> - (A-A„)V55Agt ' ^'^^^''^O' - -(A-Ao)2(32A^t)' 

y(Ao; t) := and cf (Aq) € 5(M>a,/;C), 1 < i < 4. As wih be shown below, ly^ 

gives rise to the leading-order {0{t~^^'^)) term: towards the proof of this statement, one pro- 
ceeds by estimating the difference, /q— /i/2- Recall, first, the following inequality, |exp{(-)} — 

1| < |(-)l sup,g[o,i] Iexp{s(-)}l; hence, |exp{A,^}-l| < |A,^| sup^^jo ^ |exp{sA,^}|, where A,^ = 
2A,(Ao + i e {1, 2}. Since, as shown in 0, 1 1 {6{-))^^ \ \^^^C;C) < ^^o™ ^^f" 

initions of /o±(Ao), p±{^o), and A^, i G {1,2}, it follows that sup^g[o,i] |6xp{sA^}| < 00; 
furthermore, using the Lipschitz property of ln ^j^_|[f^^^^-*jp^ , |A| < Aq, and the fact that 

r(A) G 5(r;C) and I kl l^cxj^j^.^) < 1, via an integration by parts argument, one deduces 
that 

l^b^l ^ Kl{Xo)g+K^^{Xo)glng+K^^{Xo)elnt |^b| ^ ^l(Ao)g 



with i^|'(Ao)G£°°(M>M;M>o), Similarly, one gets that, 
i£.3cxp{-2fi). ^, „3 -s(l-i)^' 



Xi(Ao)g- 



< 00 



|exp{^} - 1| < 4, -P l-P{^}l - 



< 00 



with irj(Ao) G /^°°(IR>j\,f ;K>o)) ^G{1;2}. Although the expression for the difference, /q— /i/2) 
contains many terms, estimations for the respective terms are analogous. Consider, say, the 
bound for the term corresponding to (^i, ^2, ^3, ^4) = (1,0,0, 0) in /o,a, which is denoted by 



-'0,a- 



4, :=Ya{XMt)J,"nU'^MPii'^)8'"'e-^'/Wdg. 



V2 '^'uyFiv ^ 

Using the fact that < i/< fmax := — ^ ln(l — sup^^j^ |r(A)p)<cxD, and recalling the defini- 
tions of s(Ao) and 7^(Ao), one gets that, for A G C\AA(Ao;eo), l^«(^'^o; t)| < 5^^^^^=; 

hence, letting the upper limit of integration tend to +00 (for brevity, the following notation 
is used: for matrices A and B, the inequality 1^41 < \B\ means that \Aij\ < \Bij\ 

lVl<l^l^Z^^71I5r/o Ko(^;Ao)||pi(^)|e ^1 a+de. 

In 0, it was shown that Uqq := ||j»oo( ^ ^y^"^ ''^Q^Il£^(R>o-M2(C)) hence, from this esti- 

mate and the Cauchy-Schwarz inequality for integrals, 

■(■)e~^ ■ 



Recalling the estimate for lexpjA^} — 1|, one shows that, 

l|pi(^)exp{-(.)V2}||,.(M,o;C)^^^ 
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where K^Xq) ejC°°{R>M-M>o); hence, uniformly for AeC\A/'(Ao;eo), 

I 1 I i^»(Ao)lnf 
l^0,al ^ |A-Ao|A^t ' 

— t io^ exp(— ^21) 

with i^{(Ao) G>S(]R>m;-^2(1^>o))- Similarly, recalling the estimates for |exp{ sVIaV* -"^ 

— 1|, |exp{ ^4^4^ } — 1|, and |exp{A2}— 1|, and using the triangle inequality for >C^-norms, one 

shows that the remaining terms for /o,a arc of the type 0{t~~) and 0{t~2 Int), 2<m<5, 
3<ra<6. Estimating the remaining terms of I0—I1/2 analogously, one shows that, as 
and X— 00 such that Ao>M, uniformly for AGC\AA(Ao;eo), 

where ll^2(-;^o)l!^oo(C\Ar(Ao;eo);M2(R>„)) < ^4 ^ 5(M>m;M2(M>o)), and, as A^oo, 

^«(A;Ao) -0(^(Ao)|A|-i), with fc«(Ao) ecS(M>Af;M2(K>o)). 

Repeating the whole of the above analysis mutatis mutandis for each term on the right- 
hand side of the original integral expression for x'^W which appears at the very beginning 
of the proof, one shows that, as t ^ +00 and x — > —00 such that Aq > M, uniformly for 
AGC\UKe{o,±Ao}A^(^;eo), 

\^c^X) - X5/2(A)| < (^^(A;Ao)+4(A;Ao))lnt ^ 

where h^{X;Xo) := Ei'e{o,iO}e/t (A;Ao)|r(/')|, ^2"(A;Ao) ■= T,iG{±\o}^t i^'^^o) ^ the func- 
tions e;1"(A;Ao) and e;'''(A;Ao) have the following property as A^oo, /i||''(A;Ao) + /i^(A;Ao) ~ 

and ef(Ao)G<S(M>M;M2(M>o)), /€{±Ao}; moreover, I (sAo)! |£oc(C\Ar{0;eo);M2(M>o)) < 
ll^^(-;^o)||^.(C\u.,,,,„jAr(H;.o);M.(M>o)) < ^ti^r) e /:-(M>m;M2(M>o)), 4(A;.) g 
5(M>M;-A<f2(IR>o))) and Xi/2i^) represents the sum over all /1/2-hkc terms in which the 
upper limits of integration tend to +00. One can write Xi/2(A) the following form. 



where 



n,(\r.\- 1 f=^Ao. /l-|r(.?)p\ d<; Ti , ( - ln(l-r(.;)K?)) d<; 



'^11 l.'^/ — 2^"^ Jo (A-Ao) 27rj Qp^ JO (A-Aq) 27ri 

i- 2^'" Jo (A+Ao) 27ri 2^^^^^ Jo (A+Ao) 2m ' 

^-Swx^ - ^b(Ao) rsT ^n^W^^'-^^'"' , rH(Ao)e^- rej ^^^iW^f^l! ^ 

^12 l^J — 2-2»'' Jo (A-Ao) 27ri "I 2^^^^^^ Jo (A-Aq) 2-Ki 

2-^*" Jo (A+Ao) 2m ¥^ Jo (A+Aq) 2m ' 

'^21 l.'^; — 22"- Jo (A-Ao) 27ri Jo (A-Ao) 2m 



-r -^luT- Jo (A+Ao) 27ri ^2^5:^ Jo (A+Ao) 

y^r(^^ _ rB(Ao) m <2l'^W^'''^^'"' , rB(Ao)e2-- <f i" Wg^'-^^"'^' d, 
^22 K"^) — 2-^*^ Jo (A-Ao) 27ri F^^^^ Jo (A-Ao) 2m 



2m ' 
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2-2»'^ JO (A+Ao) 27ri "I F^^^^^ JO (A+Aq) 27rj ' 



rB(Ao) := r(Ao) 11 

«=n.+l 



(Ao-A;)(Ao+A;) ' 

(A„-AO(Ao+AOj ' 



■ 2 -2 
. 2 



' 2 f 



P2I 

■ 2 



' 2 Te T' ' 2 Te^' 

A-f^lf (,) = A-gno, A-^I^CO = -^I^^i^), 

'^±,21^ = ~'^±,2r('^)' '^±,22 ('^) = -^±,22^ ('>")' 
Pl2 — \Tf ■ \ ' /^21 — P12, 

the integrals are evaluated along the rays (0,efc) (and their complex conjugates), ei := 
ooexp{^}, £2 :=oo exp{ — ^}, D±i,y{-) is the parabolic-cylinder function |jl^, and 

A2=exp{ E (/9f(0)+pf(0))}, 
/e{±} 

2^ 

■2m' 



4(0) = -2^/o^'"lnk|dln(l-|r(,)P), /^^(O) =/°^ 

j,2 j,4 

'^u' (A) = 4'^ (0) /o" ^-^y^' 1^-4"^ (0) /o^ ^-^A^"'^' |& , 

^f2^(A) = -4+)(0) /- ^i^l- + ^W(O) Jo- 

X-c (A) = ^(-) (0) /f - (0) /f 4 , 

A 27rj ' 



A'I^(A) = -^(+^(0) j^^ -'-^f ^ 1^ + ^^+^(0) 



^(+)(0) = (-1)^-Ar{ (.^^ - (xf^}, ^r^(O) = ^[:^^(0), 

N N 
rc(0) := r(0)exp{4i ll): rc(iO) := r(iO) exp{4i ll], 

l=n+l l=n+l 

where r(0) := (r(A)|^g-[^)|A=0) and r(fO) := (r(A)|^^.j^)|A=o: the explicit expressions for 
^ ^ {1; 2, 3, 4}, J, /cG {1, 2}, are not written down here since they will not actually 
be needed. Since § (r(A)|^^j^)|A=o = (KA)|;,g,M)U=o = 0, ^J:^^ (0) = ^^"^(0) = /i+(A;Ao) = 0; 
hence, xj^-'^ (A) = 0, i, j E {1, 2}. To obtain the expression for Xj^-"^ (A), i, j € {1, 2}, given above, 

use was made of the explicit representation for /xoo(-;Ao) '■= (Id — C'oo(-;Ao))~^I on L'uL' 
(recall the definition of I1/2) in terms of parabolic-cylinder functions given in Sec. 7 of 

Now, substituting the expressions given above for '^±_fj'^(?) and Xj^^^^{q), jG{l,2}, into 

the corresponding integrals for Xj^''{\), i,jG{l,2}, and using the following identities p9t| , 

d,Da{<.) = \{aDa-i{^)-Da+i{^)), (?) = ^+1 (O+^i^a-i (?) , and |r(ii/)|2 = 7r/(z.sinh^z.). 
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as well as the following integral [|T9| , 



/o°°exp(-^)x"-il)_,(x)dx 



v^exp{- i (a+b) In 2}T{a) 



10 "^^i^y 4 
one obtains the result stated in the Lemma. 



r(i(a+fe)+i) 



5R(a) > 0, 



Proposition 4.2 A ^ oo, A G C\ (f U {Uf^^{{±Xi}U{±Xi}))), x(A) has the following 
asymptotic expansion, 



X(A) 



+ |q>^(x, t) + 4|: - ) ) + o(A-2^ 



I sa 

raic equations, 



(94) 



where lim t;A))i2 :=(5'^(x, t)/2A, {aj, Wjjf^]^ satisfy the following non- degenerate sys- 

A— >oo 

iem o/ 2n linear inhomogeneous 








" 9i+x52(Ai) " 






5^X12(^2) 


a™ 




5nXl2(An) 






5^Xfl(Al) 


CJ2 




5^X11 (A2) 


_ t^n _ 




ffnXll(An) 



(95) 



where, for i,j G {1, 2, . 
follows. 



M - 



■ ,n}, the nxn matrix blocks, A'^ , 13^, , and , are defined as 



+g+ Xf 2 ( )xi2 ) 3.+ W(x; 2 ( Ai ) ,xi2 ( A, ) ) 

2g+ i~^^X22 i^i)xl2 i^J )+>^J X^2 ( Aj )X^2 ( )) 

Xi2(A,){A^A^) ' 
_ 2g+ (Aixi2 (Ai)xf Alxii (A7)xf2 (Ai )) 

xJi(a7Ka?3a7') _ ' 

_ 25+ ( Aixi2 ( Ai )X^ i(_Ai )- A,_x^i ( Ai )x? 2 (Ai )) 
X?i{A,){A2-a/) ' 



« = J, 

^ / 3, 
i = j, 



CI 



29+(-AiX^i(A0x^2(A_i^+A>X^2(A>)x?i(A0) 

_ _ x_^2(a.)(a7'-a2) 

23+ (- A.x^i {Ai)x?2 (A,_)_+Aj (A, )xji (A, )) 
Xi2(A.){A7'-A2) ^ 



^ = J, 
^ / J, 



A.-9+Xii(A»)x^i(A_^+A»g/-W(xii(A0,x^i(A,)) 

_ AiXii_(Ai)^ ^ ^ ■ 

_ 2g+ (Alx^i (A7)x?_i^A7)- ATx^i (A7)xfi (A7)) 

x^i(a7)(a7'-a7') ' 
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{A)^i}?=i satisfy the following non- degenerate system o/2n linear inhomogeneous algebraic 
equations, 




' Pi ' 










92X22(^2) 












5^X21(^1) 






52"X2l(^2) 









(96) 



where, for i,jE: {1, 2, . . . ,n}, the nxn matrix blocks, S'^ , J^'^ , Q'^ , and H'^ , are defined as 
follows, 

^1X22(^1) 



£^ - 



29+ (A»xj2 (A»)xii(Al)- Alx^i (^)X^2 

x!i(a7Ka2-a-') _ ' 

29+ (AixJa )x^i^^, ) - A,_Xn ( A, )x^2 (^i)) 
X?i(Ai)(A2-A/) ^ 

2^(AiX?2(^»)x^i(A7)-A7x^i(Al)xE2(^i)) 

_ x^2(^«K^'^^?) _ ' 

29+ (A,X^2 (A,- )xii ( A. K A, xj 1 ( Ai )xi2 (Ai )) 
Xi2(A.)(Ai'-A2) ' 



^ 7^ J, 



with 



A»+g+x^i (Ai)x^i (^_»)-^» gj W(x^i (AQ.xii (A,)) 

A,Xn(Ai) 

2g+ (a7x?i (^)xii (a7)- A7x^i (A7)x^i (AI)) 



X^i(Aj){Ai'-A/) 



« = J, 



gt := Cje 



2i\'^.x+iiXH 



N 



(Aj-A;)(A,+A;)y 



1 < i < n, 



5^ {Xk;X{j) , A; G {1, 2, . . . , n}, (/iwen in Theorem 2.1, Eq. (18), and W{xi^{z),Xi'j'{z)) is the 



Wronskian of xtjW and x-,,-,(A) evaluated at z {i,j,i',j' G {1,2}); W(Xi,(z), Xi/,/(2)) : 



{xU>^)dxxh^{X)-Xh>WdxxUmx-- 
Proof. For l<z<n, set 



res(x(A);Ai 



CX^ Oj 



res(x(A);A,)= ? | 



(97) 



Prom Eqs. (92), (93), and (97), and the polar (residue) conditions in Lemma 3.3, one gets a 
system of linear algebraic equations for {cKj, Oj, bi, Ci, di, Wj, 6i}f^i: from this system, one 
shows that, for l<i<n, 



(Xl2(^i)"i + X22{^i)0'i)9t = ^ Oi 

(Xi2(Ai)A + X22(A0^*k+ = ^ 6i 
(Xii(Ai)ci + X2i(AiM)5i" = ^ Cj 
(xfi(A^)di + X2i(A^)^i)5? = ^ 



.xi2(Al) 

xCl^ 



.X52{A 



'xi^iAiy 

.XIi(Ai) 



Xfi(^ 
X^1 (A. 



Xfi(A. 



(98) 
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using Eqs. (98), which show that the matrices {res(x(A);Ai)}f^]^ and {res(x(A);Ai)}^^^, re- 
spectively, are degenerate, one simpUfies the resulting system of linear algebraic equations 
for {ai, /3i,ai,bi,Ci,di,uJi,Si}'^^^ and obtains Eqs. (95) and (96): the non-degeneracy of sys- 
tems (95) and (96) is a consequence of the unique solvability of the original RH problem 
(Lemma 2.2). Substituting Eqs. (97) into Eq. (92) and defining Q^{x,t) as in the Proposi- 
tion, one obtains, from Eq. (93), Lemma 4.1, and Eqs. (98), Eq. (94). ■ 



Corollary 4.1 As t^+oo and x— >— oo such that Ao>M and {x,t)e^n, 

Q{x, t) = Qj(a;, t) + 4i E (c^j - l^^i) + OiC{Xo) exp{-abt}) , 

where Q^{x,t):=iQ^{x,t), a and b are given in Lemma 3.3, and C(Ao) €>C°°(]R>m;C). 
Proof. Since, from Lemma 3.2, Q{x,t) = 2i lim (Am''(a;, t;A))i2, the result follows from 

A— >oo 

Lemma 3.3 and Proposition 4.2. ■ 
Proposition 4.3 Ast^+oo and X — ^ — DO such that Xq ^ M a/id (x, t) G Jl^j, 

Qix,t) = Q<iix,t) + 4i{un - |^«n) + + 0{C2iXo)e-^°'°'), (99) 

where 

a = si Xl^L ( An ) +^22 ffn Xf 2 ) (100) 

(aiia22-ai2a2i) ' ^ ^ 

ailffn Xn (An)+a2ign Xf2 (An) /im\ 

LOfi — 7^ — ^ ^ — ^ — ^ , liUil 

(011022-012021) ^ ^ 

;^ An+gnX^2(An)xi2(An)+Ang;tW(Xi2(An),Xi2(An)) nnO^ 

— AnX^2(An) ' ^^"^^ 

._ (AwX22 (Aw)Xll (An)~AnX21 (An)Xl2 (Aw)) (103) 

Xll(An)(A^— An ) 

;^ 2gn(AnX22(An)Xn(An)-AnXil(An)Xl2(An)) finA\ 

X22(An)(An -A;^) 

-> A„-g;tx2i(An)xn(An)+An gn W(xii(An),Xn(An)) nnc;^ 

- ^^^) ' 

Q+ix,t) is givenin Theorem 2.1, Eqs. (20)-(22) and (24), ao :=min(a, 8min{77i}^/) (> 0), 
6o:=min(6,min{|en-C;|}r=i'). Ci(Ao) G<S(M>m;C), and C2(Ao) g£°°(R>m;C). 

Proof. Solving Eqs. (95) as t ^ +oo and x^ —oo such that Aq > M and (x, t) G n„ for 
{ctiliLi and {LOi}^^i via Cramer's rule, one shows that 

ai, uji ^ 0{exp{-a^ min \U-^l\t}), 1 < i < n - 1, (106) 

l<;<n— 1 

where a!' : = 8 min{?7i}J!~|^'^ (> 0), and a„ and a;„ arc given by Eqs. (100) and (101): the result 
now follows from Corollary 4.1 and the estimates in Eq. (106). ■ 

Proposition 4.4 As t^+oo and X — > — oo such that Xq^ AI and (x,t)Grijj, 

Q(x,t)=g+(x,t) + 0(^i^), 
where Q'^{x,t) is given in Theorem 2.1, Eqs. (14)-(29), and C(Ao) g5(M>a^;C). 

Proof. Substitute Xij{')^ ■^)j€{l,2}, from Lemma 4.1 into Eqs. (99)-(105) and neglect 
exponentially small terms. ■ 



5 Asymptotic Evaluation of ((^ ^(x, t; 0))ii)^ 

In this section, the phase integral, ((^'^^(x, t; 0))ii)^, which appears in the gauge trans- 
formation (Proposition 2.3, Eq. (9)) is evaluated asymptotically as >+oo (x/t~0(l)). 



26 



Lemma 5.1 As t^+oo and X — > — oo such that Xq^AI ciTid (x,t)Gf^'n,, 



((^-i(x,t;0))n)^ = exp{21n(xi2(0))}exp{f(/, 



2ifrXo ln(l-|r(g)p) 



N 

E 

l=n+l 



g 

n 



dQ - Jo 



oo lii(l+|r(ig)p) 



dQ)] 



X exp{-4i E ^;}exp{21n(l- E(f" + f^))} + 0(C(Ao)e-"^*^ 



i=l 



where bj 



and C{Xo)e£^ 



'X.22 



l<«<?i-, {/3«,^i}"=i satisfy Eqs. (96), a andb are given in Lemma 3.3, 



;C). 



Proof. From Lemma 2.1, the proof of Lemma 3.1, Proposition 3.2 (Eqs. (85)), Lemma 
3.2 (Eqs. (86)), and Lemma 3.3, one gets that 

N 

*(x,t;0) = x(0)(5(0))-3 n (di+(0))-3 + 0(Ci(Ao)exp{-a5t}), 

l=n+l 

where Ci(Ao) G >C°°(M>Af ;M2(C)). From Propositions 2.1-2.3, the parametrization for the 
discrete eigenvalues (Sec. 2), Lemma 3.2 (Eqs. (87)), Proposition 4.1 (Eqs. (92) and (93)), 
the proof of Proposition 4.2 (Eqs. (97) and (98)), and the cti and £73 symmetry reductions 
for x(A), one shows that. 



^-l(x,f;0) = (/i(0))'^ 



1-E 



V 



26^ _^ 25j,- 







'2ai I 2ci 



\ 



+ 0(C2(Ao)e- 



-abt\ 



i=l 



where h{0) := X22(0)(^(0)) ^ exp{-2zX;/In+i7/}i bi, l<i<n, and {Pi, ^ are as given 

in the Lemma, q = — J- uJi, 1 < i < n, {ai,uji}f^i are defined by system (95), and 

C2(Ao) £>C°°(M>m;^2(C)): using the expression for 5+(A;Ao) given in Proposition 3.1 (and 
Remark 3.1), one obtains the result stated in the Lemma. ■ 
In order to estimate (X22(C'))^7 the following proposition and lemma are necessary: 



Proposition 5.1 Define Q^[x,t) -.=21 lim (Ax'^(a^, i;A))i2. Then 

A— »oo 



2 ( TOO ln(l+|r(te)p) 



dQ- 



00 ln(l-|r(e)|^) 




dg), 



ixhm' = is^ioMrew{^iu\QKQ,trdQ}. 



Proof. Follows from the definition of X^(A) given in Proposition 4.1, Proposition 2.2, 
and Proposition 8.1 in [[7|. ■ 



C(Ao)(lnt)2 
^372 



Lemma 5.2 (|^2[) As and such that Xq>M, 

+ /\ \ i(4Alt-i^(An) Int) + N 



where 



"ti,o(^o) = y^exp{i0+(Ao)}, 



-(Ao) = 0+(Ao)-f +argr(ii/(Ao))+argr(Ao)-3i/(Ao)ln2+2 E arg 



N 

E 

«=n+l 



(Ao-A;)(Ao+A;) ^ 

(Ao-Ai)(Ao+AO;' 



-l,2,o('^0 



'SirAg 



dg 



lexp|i(4 E 

=0/ [ V /=n+l 



7z + 2i?+(Ao) 



(/>"'"(•) zs (7i?;en in Theorem 2.1, Eq. (22), and C(Ao) G 5(IR>m;C) 



rAp ln(l-|r(g)p) dg , roo ln(l+|r(ig)p) 
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Comment to Proof. Up to the leading term, the asymptotic expansion was 

proved in Q. The 0{t~^) term constitutes the leading-order contribution from the first- 
order stationary phase point at A = 0: the complete proof of this asymptotic expansion can 
be found in [^]. ■ 

Proposition 5.2 As t^+oo and x^—oo such that Xq>M, 

{XIM? = exp|^('y||^"^/^(i?+(0)cos(K+(/.;t))-i?+(0)sin(K+(A.;t)))^]| 



+ 0\ 



t JXo /x^ 

'C7(Ao)(lnt)2^ 



Xlt 



Where i?+(0) = 9{i?+(0)}, ii+(0) = 5R{i?+(0)}, i?+(0):= ^ 



• exp{4iEiln+i 7/}; '^+(Ao; t) :=4A^t-z.(Ao) lnt+0+(Ao)-2i?+(Ao), and C(Ao) GcS(M>a/;C). 

Proo/. Writing J^JQHQ,t)\^d0 = -{\\Qi{;t)\\c2iR.C))^ + I-JQHQ,t)\^d0, using the 
expressions for *)II£2(]^.q)^ ^iid (X22('-'))^ given in Proposition 5.1, the asymp- 

totic expansion for Q^{x,t) given in Lemma 5.2, the following inequalities, |exp{(-)} — 1| < 
l(-)l supsg[o,i] |exp{s(-)}| and 0<i/(Ao)<fmax = -^ln(l-sup^g]^ |r(A)p)<oo, and the fact 
that r(A)G5(r;C), one obtains the result stated in the Proposition. ■ 
Lemma 5.3 As t— >+oo and X — > — 00 such that Xq^AI ciTid (x,t)Gf^'n,; 

ii^-Hx, t; 0))n)2 = exp{i arg q+{x, t)} + , 

where arg(7+(x,i) is given in Theorem 2.2, Eqs. (48)-(51), and C(Ao) g5(IR>a//;C). 

Proof. According to Lemma 5.1, in order to evaluate t; 0))ii)'^, estimates for 

exp{2 ln(x22(0))} ^iid required: the estimation for exp{21n(x22(0))} is given 

in Proposition 5.2; hence, it remains to estimate {bj, Si}"^^]^. Solving system (96) as t— >+oo 
and X — > — oo such that Aq > M and (x, t) G 0„ for {/3j}"^;^ and {5i}'^^i via Cramer's rule, 
one shows that, j3i, (5i ~ 0(exp{— a'' min |Cn— l<i<n— 1, and 



fin 



l<i<n-l 



(c+ ■)-/+ _ 77+ c+ 1 ' " (c+ ■)-/+ _-p+ n+ \' 



where 



q+ .^ gtxlii^n) |g;!:|"xii(An)xi2(^n) |g;tPxi2(-^n)v^(xji(A„),xii(An)) 

• X^l(An) ^ A„ X?l(An) 

2A;ri(,+|2x^2(A^)xii(A;r) 

(A2-A„") ' 

?+ -p+ _ ^+ (5^+ 1 I gn^(xil(An),Xn(An)) , gn W(x!2 (An) ,Xm (An )) 

C-nn^nn -^nn^nn • ;^c^(A„)x^i{A„) ^ X^2 (An) ^ X?i (An)xi2(^n) 

, gx^i (A;r) _ 3;tx'f2(^) , (2A;r)2|<;+px^2(A„)x^i(A;r) 

A„xi2(^n) A„X?l(An) (A2-A„^)2 ' 

?+ gnX2l(An) _ 2An Iffn PXn (An )xi2 (An) 

'^^ - X^2(An) ^71^^^^ ■ 

Substituting the expressions for Xiji')-, hi G {1,2}, given in Lemma 4.1 into the above 
equations for 5+, and £nn'^nn~-^nnSnn} s-J^d recalling that (Eqs. (98)) bi = — ^c^(a') Ai 

l<i<n, one obtains, as a result of Lemma 5.1, keeping only 0{1) and 0{t~^/'^) terms, the 
result stated in the Lemma. ■ 

Corollary 5.1 Ast^+oo and x^—oo such that Xq>M and {x , t) £ iln , 

q{x, t) = Q+ (x, t) exp{i arg g+ (x, t)} + O (^MlLll! ~ 
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where Q^{x,t) is given in Theorem 2.1, Eqs. (14)-(29), axgq^{x,t) is given in Theo- 
rem 2.2, Eqs. (48)-(51), anrf C(Ao) g5(]R>m;C). 

Proof. Consequence of Proposition 2.3 and Lemma 5.3. ■ 

Corollary 5.2 As +oo and x +oo such that Xq := — > M, f > ^, s € ]R>o, 
and {x, t) G Cl^, 

u{x,t)=v+{x,t)wUx,t) + o(^^^^), 



where v^{x,t) and vj^{x,t) are given in Theorem 2.3, Eqs. (58)-(70), and C(Ao) g5(M>m; 
C). 

Proof. Consequence of Proposition 2.4 and Corollary 5.1. ■ 

6 Asymptotics as ^— ^— oo 

In this section, the asymptotic paradigm presented in Sees. 3-5 is reworked for the case when 
t-^ — oc: since the proofs of all obtained results are analogous, they will be omitted. This 
section is divided into three parts: (1) in Subsection 6.1, extended and model RH problems 
are formulated as t^—oo; (2) in Subsection 6.2, the model RH problem formulated in (1) 
above is solved asymptotically as t^—oo for the Schwartz class of non-reflectionless generic 
potentials; and (3) in Subsection 6.3, the phase integral, {{^~^{x,t;G))\\f' , is evaluated 
asymptotically as t^—co. 

6.1 Extended and Model RH Problems 

Proposition 6.1.1 In the solitonless sector (.2^ = 0), as t— oo and x^+oo such that 
Xo>M, 

m{x,t;X) = A{X) + 0{^), 

where A{X) := {S- {X;Xo)y^ , 

<5-(A;Ao) = ((A-Ao)(A + Ao)r^exp{ E PiW}, 

i€{±} 

p±w = -h itZ - - 1^(^)1')' 

u := iy{Xo) is given by Eq. (21), ll('5"(-;Ao))=^^||£oo(C;C) •= ^up^^^ |((5-(A;Ao))=^^| < oo, 

(^(5~ (ibA;Ao)^ =(5~(A;Ao), the principal branch of the logarithmic function is taken, ln(|Lt— 
A):=ln|/x-A|+iarg(/x-A), arg(/x-A) G (-vr, tt), and C(Ao) g5(M>m;M2(C)). 

Remark 6.1.1 Hereafter, all explicit x,t dependences are suppressed, except where abso- 
lutely necessary, and 5~ {X;Xq):=6{X). 

Lemma 6.1.1 There exists a unique solution m^(A) := m{X){A{X))-^:C\ {Z^ U T) ^ 
SL(2,C) of the following RH problem, 

1. m^(A) is meromorphic VAGC\r, 

2. 

fh^(X) = fh^{X)v^{X), A G f , 

where 

^A^^ ^ ,-^e(A)ad(a3) I -_r(A)Kl))5- (A))-i KA)L(A)5+(A) \ 
V -r(A)(5_(A))-i(5+(A))-i {5.{X))-H^{X) )' 
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3. m (A) has simple poles at {ztAj, ibAi}^^ with {l<i<N) 

res(m^(A);Ai) = lim m^(A)i;i(5(Ai))-V_, 

A — >\i 

res(m'^(A);-Aj) = -(T3rcs(m^(A);Ai)cT3, 
res(m^(A);A^) = \im_m^{\)vl{5(\i)fu+, 

res(m'^(A);-Aj) = -cT3res(m^(A);Ai )(73, 

4. as 00, \^C\{Zdyjf), 

moreover, Q{x,t) = 2i lim (Am'^(x, t;A))i2 is equal to Q{x,t) in Lemma 2.2, Eq. (11). 



Proposition 6.1.2 Introduce arbitrarily small, clockwise- and counter-clockwise-oriented, 
mutually disjoint (and disjoint with respect to T) circles and , l<j<ra — 1, around 
the eigenvalues {±Aj}"~| and {±Aj}"~]^, respectively, and define 

n—l 



m\\) := { 



m^A), AGC\(ru(U(i^f UL±))), 

i=l 




vi{5{±\i))- 



(AtA,) ^-j' XemtKf, l<i<n-l, 
A G intZf , 1 < ?; < n - 1. 



Thenm\X) solves^a RH problem on (o-i;\U^r/({±Aj}U{±A^}))U(U^ri^(^^uZf )) with the 
same jumps as m^{\) on T, m^(A) = m^(A)i;^(A), and 



mL(A)(l+M^(±M^ 



IA?A-^^- 
(ATAi) + 



AG iff, l<z<n-l, 
AeZf, l<i<n-l. 



Remark 6.1.2 The superscripts it on {K^}^^^ and {Lf}i'=i, which are related with 
{±Xi}^~i and {±Xi}^~^, respectively, should not be confused with the subscripts ± ap- 
pearing in the various RH problems in this and the next subsection, namely, m±(A), m^(A), 
m^(A), mi(A), x±(A), ^±(A), and x^(A). 



Lemma 6.1.2 Set 



, n—l ^ n—l — , ^ . 

m^(A) n (c^/+(A))-'^^ \eC\{Tu{[j{Kfu Lf))), 

1=1 i=l 

n-1 ^ 

^ (A)(J;?±(A))-i n (A))-"^ AGintXf, l<i<n-l, 

1=1 

~ n—l ^ 

m^(A)(J7±(A))-i n (^^^_(A))-'^^ A G intL±, 1 < i < n - 1, 



where 



dl,{\Y= %Z'i!^^X% AGC\(UJiff ULf)), l</<n-l. 



n-l 



± I I r ±^ 



a!i_(A):=<^ 



AGUintKf, l</<n-l, 

(a±a7) 

(A-A,)(A+A,) 



i=l 
n-l 



, A G U intLf , 1 < / < n - 1, 
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and the SL(2,C)-valued, holomorphic in intK^ and intL^, respectively, functions J^±(A) 

i 

and J^±(A), l<i<n—l, are given by 



-i;i((5(±Ai)) n 



V '71 

/ n-1 dr^A) 



C« n-1 d-_i(A) 



(A) ^(«(±a7))2c« (A) 



n 



rf/ I (A) (rf, (A)) 



(A) 



(A=FAi) 



with 



n-1 



C" = (^i)-'(<^(±A0)2(di_(±Ai))2 n (ciz+(±A0)2, 1 < i < n - 1. 



i=i 



T/ien m«(A):C\((Zd\U^r/({±Ai}U{±Ai}))U(ru(U7r/(Kf UL±))))^SL solves the 
following, extended RH problem on (ct^ \ U"r]^({±Ai}U{±A^}))U(U"r/(Kf ULf )), 

ml{X) = m"_(A)e-^^W^'i('^3)??tt(A), 

where 

I — ~ n-1 



\ 5-{A)5+(A) \^(A-AO(A+AO 



r(A) 



(5_(A)5+(A))-i UA-AO(A+Az) 



(A-A;)(A+A;A^ \ 



,3+ (A) 
5_(A) 



?«(A): 



^+ (^) i^A— j n I ApAf ) A G U , 



2Ai 



j=l 

n-1 



i=l 



mi/t po/ar residue ) conditions, 



^ n—l / \ 2 

res(m«(A);Ai) = \\m fh^{\)vi{6{\i))-^ ' i^i^MMM 



o"_ , n < i < N, 



res(m''(A);— Aj) = — (73res(m''(A);Aj)(T3, n < i < N, 

^ n— 1 / — — — — \ 2 

res(m«(A);Al) = lim m«(A)w(5(Al))2 n ?=^rWTl4 " < ^ < ^> 

A^Ai i=l V^^^'-^'^^^'+^'V 

res(mS(A);— Aj ) = — (T3res(m''(A);Aj )(T3, n < i < N, 
and, as X^oo, AGC\((Zd\U,"rii({±Ai}U{±Al}))U(ru(U,t"/(i?^ uZf)))), 



fl(A) = I + 0(A-i); 



moreover, Q{x,t) = 2i lim (Am^(x, i;A))i2 is equal to Q{x,t) in Lemma 2.2, Eq. (11). 

A— >oo 



Lemma 6.1.3 Let x(A) solve the following RH problem on cr£\[J^^-^ {{±Xi}U{±Xi}), 



X+(A) = x-(A)e 



-i6l(A)ad((T3) 
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with polar ( residue ) conditions, 

res(x(A);A,) = Y^mx{\)vm\))-^ "tf (g';:^jg";t^j) > n < ^ < iV, 
res(x(A);-Ai) = -(73res(x(A);Aj)(73, n<i<N, 

n— 1 / — — — — \ ^ 

res(x(A);Al) = ^lim x(A)w(5(Al))2 U [^E^^j^) <^+, n<i<N, 

res(x(A);-Ai ) = -cr3res(x(A);Ai )a3, n<i< N, 
and, as X^oo, AeC\(U,Cn({±Aj}U{±A7})uf ), 

X(A) = I + 0(A-1). 

Then as t ^ — oo and x — *■ +00 such that Aq > M and (x, t) G Jl^, i/ie function E{\) := 
m''(A)(x(A))~''^ has the following asymptotics, 

E{X) = I + 0(F(A;Ao) exp{56i}), 

where \\nM\\c^^QM.iC)) < ^' •)IU(M^^;m.(C)) < ~' F(A;Ao) ~ o(^) as 

A^oo TOf/i C(Ao)e>C°°(R>M;M2(C)), 5:=8min{r?i}^ri^ (>0), and 6:= min{|^„-,ei I }r=i^- 

6.2 Asymptotic Solution for x(A) 

Proposition 6.2.1 T/ie solution of the model RH problem formulated in Lemma 6.1.3, 
x(A):C\(fu(Uj^„({±Ai}U{±A^})))-^SL(2,C), has the following representation, 

~ [ X-(g)(^"(g)|f-I) dQ 
X(A) = X.(A) + — , (107) 

where 

~, tW —T \ ^ ( res(x(A);Ai) o-3res(x(A);Ai)CT3 , res(x(A);Ai) CT3res(xW|Ai]CT3 \ 
Xd{^}-i- + ^^J^ (A-AO (A+A,) + (A-A-) (A+A7) j' 

The solution of Eq. (107) can he written as the following ordered product, 

x(A)=Xd(A)x^(A), 

where Xd{X) is given by Eq. (108), and x'^(A) solves the following RH problem: (1) x'^(A) 
is piecewise holomorphic VAGC\r; (2) x+(A) =X-(A) exp{— i^(A)ad(c73)}(t;''(A)|j^), AgT; 
and (3) as X^oo, AGC\r, x^'lA) = I+0(A-^). 



Lemma 6.2.1 Let ?o denote an arbitrarily fixed, sufficiently small positive real number. 
For K G {0, ±Ao}, sei jV'(K;eo) := {A; |A— K|<eo}. Then as t^— 00 and x^ +00 such that 
Xo>M and A G C\U~^^Q _|_^^-^A/"(H;?o), x'^(A) has the following asymptotic expansion, 

+ exp{i(0- (Ao) + (Ao; t))}a+) + O 

where z^(Ao), ^~(Ao), and $~(Ao;t) are given in Theorem 2.1, Eqs. (21), (23), and (24), 
l|G(-;Ao)IU(C\u~ Ar(K;7o);M.(C))<~' G(A,)g5(R>m;M2(C)), G(A;Ao)~o(^) 
as A— ^00 with C(Ao) G (S(]R>m;-^2(C)), and satisfies the following involutions, x'^(— A) = 
o-3X''(A)o-3 and x''(A)=crix'=(A)cri. 



G(A:Ao)ln|^| 
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Proposition 6.2.2 For n<i<N , set 
res(x(A);Ai) = 



Pi K 



res(x(A);Ai 



d7 S7 



Then as A oo, AG C\ (FU (U^„({ibAj}U {itAi}))), x(A) has the following asymptotic 
expansion, 



TV 



X22(Ai) 



a,- 



where lim t;A))i2 := Q^{x,t)/2X, {a- ,lo^ }iLn satisfy the following non- degenerate 

system of 2{N —n+1) linear inhomogeneous algebraic equations, 

_5nXl2(An) 
5'n+lXl2(An+l) 

5WXi2(Aiv) 

5nXii(A„) 










«n+l 





















5n+lXll(An+l) 

where, for i,j e {n,n+l, . . . , N}, the {N—n+l)x{N — n+l) matrix blocks, A", B~ , C" , 
and T>~ , are defined as follows, 

^1X22 (-^«) _ „ 

2g,- )x^2 (Aj )+ Aj x^2 (Aj )x^2 i^i)) 

X^2(A,)(A2-A2) 



B. 



(A,X^2 (>^^)x•i-,_<^^ ) - AiXil ( Aj )xf2 (^i )) 

_ 25- (AixS2 (A«)xJi(Ai )- AiXii (A,)xf2(Ai)) 
xJi(A,)(A2-a/) ' 

_2^(AiX-2(A.)x?i(A7)jA7x^i(A7)xj2(Ai)) 

_ x_^2(Ai)(A,"-A|) _' 

2gr (- AiX^i (Ai)x^2 (A^+A, Xi2 (A, )x^i (Aj )) 
Xi2{Aj){Ai -A|) 

Al-^X^i (ApX^i (Al)+Al^W(x^i (A7),xfi (M)) 

_^»^ii(^'L^ ' 

_ 2flr (Aix^i (A.)x?i^A,0- A.^^i (A, )xJi (Ai)) 
X^i(A,)(A.'-A/) ' 



« = J, 

i = j, 
i = j, 
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{Pi , }iLn satisfy the following non- degenerate system of2{N—n+l) linear inhomogeneous 
algebraic equations, 










Hn+l 




Q 1 1 Vool Anj_i ) 












5nX2l(^n) 






5n+lX2l(^n+l) 


- 




9nX2i(^n) 



where, for i, j G {n,n+l, . . . , N}, the {N—n+l)x{N — n+l) matrix blocks, £ , !F , Q , 
andTi", are defined as follows, 

( Ai-g-xf2(-^i)x^2(^»)+^»grW(x^2(^»).Xi2(Ai)) 

-^iX22(-^') 



* = J, 



2g,: (^jX^2 (WX22(^i)-^iXi2(>'i )X22 i^j )) 

Xi2(A.){A?-A2) 

x|i(a7Ka2-a7 ) 

2g," (AiXf2 (AQxiiXAj )- A^^i (Aj )x§2 (Aj)) 
Xfi(A,)(A2-A/) 



, i = 3, 



25r(AiXf2(A«)xii(Ai)-Aixfi(Ai)xi2(Ai)) ^ ^ • 
_ X^2(Ai)(A-'-Af) ' 

2gr(AjXf2(Aj)xii(Ai)-A,xfi(Ai)x§2(Aj)) 



Xi2(Aj)(Ai -A|) 



Ai+g-x^i(Ai)xii(A_i)-AigrW(x^i(Ai),xSi(Ai)) 

AixJi(Ai) 

2g- (Alxfi (ADxii (A7)- A7x^i (A7)x^i (AI)) 



x^i(a,){a7'-a/) 



(7~ =-^12(Ai) - A- =_X^2(Ai)^- - 

Xi2(Ai)"' ' Xi^CAi)^^ ' 



_ X^i (A») 
X?i(Ai) 



« = J, 

_ X^i (Al) 



< i < AT, 



5 (Afc;Ao), G {n, n+1, . . . , A^}, ^iwen in Theorem 2.1, Eq. (19), andW{xij{^)-,x'iiji{^)) is 
the Wronskian of XijW and Xi>.j>{^) evaluated at z {i,j,i',fe{l,2}):W{xiAz),Xi'i'{z))-= 



Corollary 6.2.1 As t^—oo and x^+oo such that Ao>M and {x,t)e^r, 



N 



Q{x,t) = Q^_ix,t) + 4it^i^u;T - ggjlaT j + 0(C(Ao) exp{a6i}), 
where Q^{x,t):=iQ^{x,t), a and b are given in Lemma 6.1.3, and C(Ao)G>C°°(]R>m;C). 



Proposition 6.2.3 As t^—oo and x^+oo such that Xo>M and {x,t)e^r: 



Q{x,t) = Q^_{x,t)+Ai(u- 



Xl2(An), 



Xi2(An) " 



Ci(Ao)ln|t| 
t 



+ 0(C2(Ao)e««^»*), 
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where ai ,cjj ~C'(exp{a'' min \U-^l\t}), n + l<i<N, a" :=8mm{r/i}/L„, ^ (>0), 

n+l<l<N 



l.''-ll'*22 '^12'*21> 



, _ Clign Xn(^n) + a2ign X^2(-^'0 
v'*ll''22~''l2''21^ 

._ 2g;:(A„X22(^")Xn{^n)-A„X2i(^n)x5:2(-^")) 
■ Xji(A^)(A2-A^') ' 



2g„ (A„X22(-''n)Xn{An)-A„X2l(An)X?2(-'^")) 

■ xi2(A„)(A^'-A2) ' 

An-gnX2l(An)Xn(An) + An ffn W(xii (An ) ,Xn ( An )) 

is ffiuen in Theorem 2.1, Eqs. (20), (21), (23) and (24), ao := min(a, a'') (> 0), 
bo : = min (6, min{ | ^„ - ) , Ci ( Aq) G 5(M>m ;C) , anrf C2 (Aq) G (M> a/ ;C) . 

Proposition 6.2.4 Ast— >— 00 and a; —> +00 sitc/i t/iai Xq>M and {x,t)Gil.n, 

Q{x,t)=Q-,ix,t) + 0{^^^), 
where Q~g{x,t) is given in Theorem 2.1, Eqs. (14)-(29), and C(Ao) €5(M>a/;C). 

6.3 Asymptotics of ((^^^(x, t; 0))ii)2 as 

Proposition 6.3.1 Define Q^{x^t):=2i lim {\x^{x,t;\))i2. Then 

A^oo 

ixhm' = r (0;Ao))2 exp{z IQHq, t^dg}. 



Lemma 6.3.1 ([^]) Ast—>— 00 and x—i- +00 such that Xq> M , 

7l 

(-t)3/2 



where 



^la.ol-^o) = y exp{i6' (Ao)}, 



g-(Ao) = r(Ao) + f-argr(iKAo))+argr(Ao)+3KAo)ln2+2garg (j;:°:j^g 



. (dir(g)\ ra) d[r{g}\ ja) \ f / n-1 \^ 



roo ln(l-|r(e)p) dg 



Mo g 

(?!)"(•) is given in Theorem 2.1, Eg. (23), and C(Ao) g5(]R>a/;C). 
Proposition 6.3.2 As t^—00 and x^ +00 such that Xo>M, 



ixUO)? = exp<j^UA|°'^/^(i^-(0)cos(^-(M;^))-i?;(0)sin(^-(^;^))^^ 



+ O 



t Jxo 
/ C(Ao)(ln|t|)2 ' 



35 



where R^{0) = 'i{R-{0)}, R-{0) = ^{R-{0)}, R-{0):-- 



' d{r(g)\ 



Q=0 



dg 



exp{4i Er=i (Ao; t) := 4A4t+z.(Ao) ln|t| +0- (Aq) -2i9- (Aq), and C(Ao) G cS(M>a,;C) . 



Lemma 6.3.2 As t^—oo and x— > +cxd such that Xq> M and {x, t) € 



((vI'-nx,t;0))n)2 = exp{21n(xi2(0))}exp{f/,^Mi^tMid^,}exp{-4i t 7^} 



n-l 

E 

X exp{21n(l- E(^ + ^^))} + 0(C^(Ao)exp{S6t}), 



^ ,2b- . 25- 



where (X22(0)) ^■^ ^iwen m Proposition 6.3.2, = > 'n < i < N , fi- ,6^ 

0(exp{a'' mill \^n-^j\t}), n+l<j<N, 



n+l<j<N 



/Q- _ Pn X- _ On 



with 

raM- .^ dnXhi^") I lgnPx2i(^)X22('^") Iffn P X22 )^ (Xn (^) 'Xgl (^)) 
■ Xfi(A„) ^ A„ x!i(An) 

2A;ri9;rPx22(An)x2i{A^) 



•Vt"(X^l(An),X^l(An)) , g-H/(x^2(An),Xi2(An)) 

^nnrinn -T^^y^^ . (x„)^c^ (x„) ^ xi2(Au)x^i(An) ^ Xfi (An)xi2 (An) 



I £ri_X2l(An) _ gnX^2(^) I (2 An )^ |3n Px22 ( An )Xn (An ) 

A„xi2(An) A„xfi(An) (A2-A;r^)2 ' 



._ ffn X2l(An) _ 2A„|g„ PXn(An)X22(An) 

" ■ xi2(An) (a^-a;:^) ' 



and C(Ao)G/:°°(ffi>M;C). 



Corollary 6.3.1 Ast^—oo and +00 such that Xq>M and {x,t)£il.n, 

q{x, t) = Q^ix, t) exp{z arg q-{x, t)} + O , 

where is given in Theorem 2.1, Eqs. (14)^(29), argg~(x,t) is given in Theo- 

rem 2.2, Eqs. (48)^(51), and C(Ao) g5(M>m;C). 



Corollary 6.3.2 As t^-00 and x^-00 such that Aq := y i(f -7) >M, j>^, seR>o, 
and (x, t) £ Qn, 

u{x, t) = v-{x, t)w^{x, t) + f ^Cxo)Mt\)A ^ 



where v^{x,t) and w^{x,t) are given in Theorem 2.3, Eqs. (58)-(70), and C(Ao) G 5(IR>m; 
C). 
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